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Abstract. We study global behavior of small solutions of the Gross-Pitaevskii 
equation in three dimensions. We prove that disturbances from the constant 
equilibrium with small, localized energy, disperse for large time, according to the 
^vq | linearized equation. Translated to the defocusing nonlinear Schrodinger equation, 

this implies asymptotic stability of all plane wave solutions for such disturbances. 
We also prove that every linearized solution with finite energy has a nonlinear 
solution which is asymptotic to it. The key ingredients are: (1) some quadratic 
transforms of the solutions, which effectively linearize the nonlinear energy space, 
• (2) a bilinear Fourier multiplier estimate, which allows irregular denominators 

due to a degenerate non- resonance property of the quadratic interactions, and 
(3) geometric investigation of the degeneracy in the Fourier space to minimize its 
influence. 
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1. Introduction 

We continue the study (13, [18] of the global dispersive nature of solutions for the 
Gross-Pitaevskii equation (GP) 

iipt + Atfj = (\ip\ 2 - l)ip, ijj:R 1+3 ^C (1.1) 

with the boundary condition lim|a;[_>.oo ip = 1. The equation itself is equivalent to 
the defocusing nonlinear Schrodinger equation (NLS) by putting ip = e~ lt ip 

icp t + Aip = \<f\ 2 (f, (1.2) 

but the nonzero boundary condition has nontrivial and often remarkable effects on 
the space-time global behavior of the solutions both in dispersive and non- dispersive 
regimes, on which there has been extensive studies [U El [31 [5j El [8j [9j [HI D21 13 [HI [T71 
HE! EH [25l E6j EH [UJ. This boundary condition, or more generally lirn| a; |_ ! . 00 \ip\ = 
1, is natural in some physical contexts such as superfluids and nonlinear optics, 
or generally in the hydrodynamic interpretation of NLS, where \ip\ 2 is the fluid 
density and the zeros correspond to vortices. (For more physical backgrounds, see 
PQ GDI E31 El E2] and the references therein.) 

Hence it is not surprising that our GP equation, after some transformations (13.151) . 
is very similar to a version of the Boussinesq equation 

u tt -2Au + A 2 u = -A(u 2 ), (1.3) 

which is a model for water waves and other fluid dynamics. More precisely, this 
equation can be regarded as a simplified version of our transformed equations (See 
Remark 13.11 for the details). Our analysis will reveal that the linear part, which 
is exactly the same for (11.31) and our GP, has better dispersive properties than 
both the low frequency limit (the wave equation) and the hight frequency limit (the 
Schrodinger equation) for this type of bilinear interaction. 

We recall some known facts about GP (11.11) . For any solution ip — 1 + u, we have 
conservation of renormalized energy: 

/ |T7 , |2 (H 2 -l) 2 , f ,2 (2Rew+ \u\ 2 ) 2 , ^ AS 

E 1 (ip):= / W\ 2 + ^-±- — J -dx= / |Vw| 2 + ^ 1 1 ; dx, (1.4) 



provided it is initially finite, and then the solution tp is unique and global [12] . As a 
remarkable feature of GP, there is a family of traveling wave solutions [31 E] of the 
form ip(t, x) = v c (x — ct) with finite energy for < |c| < v^2, where a/2 is the sound 
speed. They are quite different from the solitary waves of the focusing NLS in their 
slow (algebraic) decay at the spatial infinity [16j and bounded (subsonic) range of 
speeds [T5] . 

Recently it is proved pQ that there is a lower bound on the energy of all possible 
traveling waves for (11.11) in three dimensions: 

S := ini{Ei (ip) \ ip(t,x) = v(x — ct) solves (II. ip for some c} > 0. (1.5) 

If one believes that there is no more stable structure supported by (jl.ip . it is natural 
to expect that the regime below the threshold £ is dominated by dispersion, as 
was conjectured in [1]. Understanding the effect of dispersion in the nonlinearity 
is an essential step towards investigating asymptotic stability of traveling waves. 
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However the problem is not quite easy, because we have quadratic interactions of 
the perturbation u = — without any decaying factor, due to the nonzero constant 
background: 

iu t + Au -2Reu = u 2 + 2\u\ 2 + \u\ 2 u. (1.6) 

They are much stronger, for dispersive waves, than cubic interactions and quadratic 
ones with decaying factors, which typically arise in the stability analysis of solitary 
waves for NLS. In fact, if one considers the NLS with general quadratic terms as a 
model equation, 

iu t + Am = \ + u 2 + Ao|m| 2 + A_w 2 , (1.7) 

then the asymptotic behavior for t — > oo of solutions from small localized initial 
data is known J20J HH] only if Ao = in three dimensions. From the technical view 
point, the quadratic power in three dimensions corresponds to the so-called Strauss 
exponent, because the LP decay estimate 

||e l * A ^||x3(R3) < |^r 1/2 |bllL3/2 (M 3 ) , (1.8) 

implies that L 3 is mapped back to the dual L 3 / 2 by the quadratic nonlinearity, with 
the critical decay order for integrability. Thus one cannot get closed nonlinear 
estimates just from the LP decay, unless one starts with given final states ([18] gives 
such a result for GP). Hence we have to take account of the oscillatory property of 
the quadratic terms. Then the term |w| 2 is worse than the others, because its phase 
by the linear approximation is stationary both in space and time at zero frequency 

C 0. 

Actually, we cannot neglect the first order term 2 Re u for large time behavior, 
so we should not expect too strong an analogy with the quadratic NLS. By the 
diagonalizing transform 

u = Ul + iu 2 i-> v = ui + iUu 2 , U = a/-A(2 - A)" 1 , (1.9) 

we get the equation for v with a self-adjoint operator H = a/— A(2 — A) : 

iv t - Hv = U(3u 2 + u\ + \u\ 2 u{) + i{2uiu 2 + \u\ 2 u 2 ). (1.10) 

The equation for U~ 1 v is3 similar to the quadratic wave equation around £ — > 

u tt - Am = |Vw| 2 , (1.11) 

if we pick up the first nonlinearity U(u\) only. This equation is known [22J to blow 
up from arbitrarily small initial data in C^°. Note that the other terms U(u 2 ) + 
2i(uiU 2 ) are even worse, enhancing the low frequency through u 2 = U^ 1 Imv. The 
wave equation behaves worse in the quadratic interactions than the Schrodinger 
(and the Klein-Gordon) equation because, beside the slower LP decay of e lt ^~^, the 
waves propagating in parallel directions are strongly resonant in the bilinear forms, 
unless they have the null structure (special coefficients killing exactly the parallel 
interactions) . 

1 Here we naturally put U^ 1 on v to map the linear energy space 7J 1 (K 3 ) for v onto that for the 
wave equation, iJ 1 (R 3 ). 
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Despite of all these observations, we can prove that the solutions for GP with small 
localized initial energy disperse linearly at time infinity. Below we denote (x) = 
a/2 + \x\ 2 and (V) = y/2 — A. H 1 is the standard Sobolev space with the norm 
IMI-ff 1 = IK^^IU 2 ! an d (x)' 1 !! 1 is the weighted space with the norm || 

Theorem 1.1. There exists 5 > such that: For any u(0) G i/ 1 (]R 3 ) satisfying 

[ (x) 2 (|Rew(0,x)| 2 + \Vu(0,x)\ 2 )dx < 6 2 , (1.12) 

we have a unique global solution ip = 1 + u of ( 11. II) such that v := Keu + illlmu 
satisfies e itH v G C(R; (x^H^R 3 )) and for some v+ G (x)' 1 H 1 (M 3 ) , 

\\v(t) - < Or 1 / 2 ), H^) { e «fft,( t ) - v+ }\\ m _> 0, (1.13) 

as t — > oo. Moreover, we have Ei(tp) = \\ (V)v + \\ 2 L2 , and the correspondence v(0) i— ► 
f + defines a bi-Lipschitz map between -neighborhoods of (x) _1 i/ 1 (lR 3 ). 

Actually we have some decay t~ e also for the weighted norm in (1 1.1 3D . but we will 
not try to specify it, since it is quite small. For pointwise decay, we can derive from 
the weighted energy estimate together with the LP estimate on e~ ttH that 

IM*)||l~ ( r3) < 0{r l ), |M*)|U~ (R 3) < o(r 9 / 10 ). (1.14) 

In view of the scaling property, the optimal (weakest) weight should be (x)^ 1 ^ 2 
instead of (x) 1 . The latter is however more convenient to estimate in the Fourier 
space, where it turns into one derivative. We can see that all the traveling waves 
with finite energy have (I1.12p finite, from their asymptotic behavior [16] as \x\ — > oo. 
For that purpose, we should not use a weight stronger than (x}~ 3 ^ 2 , hence (x) -1 is 
the unique integer choice. 

By using fractional derivatives in the Fourier space, we could reduce the weight 
slightly, but probably not to the optimal rate 1/2, since our argument needs to 
sacrifice part of the decay in some bilinear estimates with Fourier singularities. 

Those singularities arise because the linearized operator e~ ttH behaves like the 
wave equation around £ — > 0, where the parallel interactions become stronger. Since 
we know that solutions of nonlinear wave equations such as (II. lip blow up, we must 
use the difference between H and \/—2A in estimating the bilinear terms, which 
gives us some degenerate non-resonance property. The singular multipliers appear 
when we integrate on the phase in the Fourier regions for such interactions. In order 
to treat them, we will derive an estimate for bilinear Fourier multipliers (Lemma 
110.11) . which allows much less regular multipliers than the standard Coifman-Meyer 
type estimates. It will also be crucial to investigate carefully the geometric structure 
of the Fourier regions for those singularities, in order to exploit their smallness in 
volume. 

On the other hand, we also have to take care of those bad quadratic terms which 
contain w 2 — U~ x \mv. The remedy was already given in [TTl [T5] . by the quadratic 
transform 

u i-f M(u) := v+ (V)~ 2 M 2 > (1.15) 
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which effectively removes those terms from the equation. We will introduce a new 
quadratic transform, which behaves better in terms of regularity for high-high in- 
teractions. It turns out that the wave-like behavior of e~ ttH becomes an advantage 
when estimating the quadratic parts of those transforms, because it gives more de- 
cay for £ — > than e ltA . Indeed that is crucial to get the scattering in the same 
topology 1 as the initial data space. 

Extending the above result to the energy space is a supercritical problem from 
the scaling view point, and so it seems beyond our current technology (it is critical 
in four dimensions, where we have the small energy scattering [17J). However, we 
can still solve the final state problem without uniqueness. Since the energy space is 
essentially nonlinear for u, we need the explicit linearizing map M(u) to state our 
scattering result. 

Theorem 1.2. For any z + £ iJ 1 (R 3 ), there exists a global solution ip — 1 + u of 
(11. ip satisfying Ex(ip) = || (V)z + ||| 2 ( K 3) and as t — ► oo 

\\M(u)-e- UH z + \\ mm ^0. (1.16) 

The proof depends on a compactness argument, so we get almost no information 
about the initial data set for the scattering. However, the above initial data result 
Theorem 11.11 implies that it contains at least those data satisfying (11.121) . 

In (I1.16p . we can replace M(u) with v if the norm is also replaced with H 1 , but the 
quadratic part is generally not even bounded in L? (see Remark 14. II for the details). 
Hence the nonlinear effect is not really vanishing for large time in the energy sense, 
but instead the map M essentially linearizes the (renormalized) energy space 

F x := {/ £ H\R 3 ) \2Ref+ \f\ 2 £ L 2 (R 3 )}, (1.17) 

where if^lR 3 ) = {/ £ L 6 (1R 3 ) | V/ £ L 2 (M 3 )} is the homogeneous Sobolev space. 
F\ was introduced in [12] with the distance 5(f, g) defined byH 

5(f,g) 2 = ||V(/ - 0)112, + \\\\f\ 2 + Vi -\g\ 2 - 2gx\\l*. (1.18) 

More precisely, the energy space in three dimensions was characterized as 

W 6 (i + f) \eeR, feFx}. (1.19) 

As for the mapping property of M on the energy space, we have the following. Let 
for any L > and k > 0, 

F,{L, k) := {/ £ F\ \ Ex{l + /) < h\ \\f\\ L , < k}, 
H\L,n):={f eH'\ \\f\\ m <L, ||/|U.<k}. { ' j 

Theorem 1.3. (1) For any f,g £ Fx, we have 

6(f,g) 2 = ||(V)(M(/) - M(g))\\h + \\M\f\ 2 ~ \g\ 2 )\\h, (L21) 
hence the map u \— > M{u) is Lipschitz continuous from Fx to H 1 . 



2 The distance in p~2] is slightly different but equivalent to the above. 
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(2) For any L > 0, there exist k > 0, k' > 2k, and a unique map R : H 1 ^, n) — > 
Fi(2L, k') satisfying M o R = id. Moreover we have 

5(R(f),R(g)) <2\\{V)(f - g)\\ Lh (1.22) 
namely the inverse is also Lipschitz in this domain. 

Remark that the smallness is required only in L 6 , not in the energy, and the L 6 
norm decays if the solution scatters. Hence we can apply the invertible maps to 
scattering solutions and their asymptotic profiles for sufficiently large time. Thus 
we can rewrite the scattering statement ( 11.161) as 

S(u(t), R(e- itH z + )) ^ (t -> oo). (1.23) 

We can also translate the above results to asymptotically plane wave solutions for 
the cubic NLS (II .2p . First, it is easy to see that for any a > and any deN, every 
weak solution ip G L^ oc (M} +d ) satisfying \ip\ = a on M. 1+d is written in the form 

(p{t, x) = ae -^ 2 +\b\ 2 )t+ibx +l c^ ^ 24 ) 

for some b G M. d and c G K.. The above solutions for (11.11) can be regarded as 
perturbations of the plane wave ip = e~ lt {l + u), and perturbations of the other 
plane waves are generated by the invariance of NLS under the scaling and Galilean 
transforms. More explicitly, for any solution ip = 1 + u of (II. ip . for any a > 0, 
b G M 3 and c G K, we have a solution ip for (II. 2p of the form 

<p(t, X) = ae ^+\b\ 2 )t+ibx +l c^ + u ( a 2^ ^ _ 2bt ^y (L25 ) 

In particular, Theorem 11.11 implies the asymptotic stability of all plane waves (ll.24p 
in the weighted energy norm (11 . 13)) . as well as in L°° (I1.14p . for perturbations 
u(0) satisfying (I1.12p . Similarly, Theorem 11.21 implies that there exist plenty of 
asymptotically plane wave solutions with finite (renormalized) energy, at least one 
for each asymptotic profile from H l , a > 0, b G M 3 and c6R. 

Notice that the situation is quite different for the standard boundary condition 
a = 0. Indeed the invariant scaling ip i— > Xip(X 2 t, Ax) also changes the energy 

E (ip)= [ iV^f + ^Lfe, (1.26) 

homogeneously, so the smallness of the energy Eo(ip) has no impact on the scattering, 
and the transform M is reduced to the identity. Moreover, it is easy to observe that 
the free propagator e ttA is continuous on the energy space H 1 fl L 4 , but not globally 
bounded. Hence some free solutions from this space cannot be asymptotic to any 
nonlinear solution in the full energy. Once we restrict the solutions to L 2 (M 3 ), then 
the scattering problem was completely solved by [13], extending the earlier result 
[27] in a smaller space. But without L 2 finiteness, it is not even clear whether we 
can have a result similar to Theorem 11.21 in H 1 , or we should instead modify the set 
of asymptotic profiles. 

We conclude the introduction by reinforcing the conjecture in [1] into a precise 
statement of nonlinear scattering. Let Sq > be as in (II. 5p . 
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Conjecture 1.4. For any global solution if> G C(R;1 + F±) of ( 11. ip satisfying 
E 1 (ip) < £q, there is a unique z + G H l (J§?) satisfying E 1 (ip) = ||(V)^ + ||| 2 and 

\\M(u(t)) - e~ itH z+\\ mm -»• (*->oo). (1.27) 

Moreover, the map u(0) z + is a homeomorphism between the open balls of radius 
Sq^ 2 around in F\ and H 1 . 

The rest of this paper is organized as follows. Sections 2-3 are preparatory. 
Section 4 deals with the final data problem. Sections 5-13 are devoted to the initial 
data problem. In the appendix, we give a correction to our previous paper [18] . 

In more detail: in Section 2, we collect some notation and basic estimates used 
throughout the paper. In Section 3, we compute the normal forms including the 
previous and the new ones. In Section 4, we prove the final data result Theorem II. 21 
and also Theorem 1 1.31 In Section 5, we start with the initial data problem by setting 
up the main function space X and deriving its decay properties. In Section 6, we 
derive the weighted estimates for the normal form transform. Section 7 deals with 
the initial data. In Section 8, we estimate the Strichartz norms without weight. In 
Section 9, we treat the easier parts of the weighted estimate, namely those without 
the phase derivative, and also the quartic terms. In Section 10, we treat the main 
terms, namely the bilinear forms with the phase derivative, using a bilinear multiplier 
estimate allowing singularities, and assuming some bounds on the singular bilinear 
multipliers. Those bounds are proved in Section 11. The cubic terms with the phase 
derivative are estimated in Section 12. Finally we prove Theorem II. II in Section 13. 

2. Preliminaries 

In this section we establish some notation and basic setting for our problem. First, 
for function spaces, we denote the Lebesgue, the Lorentz, the Sobolev and the Besov 
spaces by L p , L P ' Q , H s ' p and Bp q respectively, for 1 < p, q < oo and s G IL The 
homogeneous versions of the latter two are denoted by H s,p and B^ q . We do not 
use the spaces with < p < 1, but instead we use the following convention 

L 1 / p = L p , L 1/p ' 1/q = L"P' q , H s ' 1/p = H s ' p } (2.1) 

for convenience of explicit computation of the Holder exponents. For any Banach 
space B, we denote the same space with the weak topology by w-B. 
We denote the Fourier transform on M. d by 



(2.2) 



ftp = : = / <p( x )e- ix tdx, 

Fl[f(x,y)] = (^ x m,y):= [ f(x,y)e**dx, 



and the Fourier multiplier of any function ip by 

rt-<v)/ := ^V(£)7(0L 

<p{-iv) x f{z,y) := {Hr l [mH[f^y)]\- 



(2.3) 



Next we introduce a Littlewood-Paley decomposition (homogeneous version) in 
the standard way. Fix a cut-off function x £ Co°(^) satisfying xi x ) = 1 f° r \ x \ — 1 
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and x( x ) = f° r \ x \ — 2- For each k G 2 Z , we define 

X k (x) :=x(\x\/k)- X (2\x\/k), (2.4) 
so that we have x k £ C£°(R d ), 

supp X k C {k/2 < \x\ < 2k}, Yl x\ x ) = 1 ( x ^°)- 



fcS2 z 

Hence the Littlewood-Paley decomposition is given by 

fc£2 z 

We denote the decomposition into the lower and higher frequency by 



(2.5) 



(2.6) 



(2.7) 



For any function -B(£i, . . . , £tv) on (M. d ) N , we associate the iV-multilinear operator 
B[fi,-.. /n] defined by 

^n.r.....fx [ B(t 1 ...fr)f 1 (t 1 )-.-fa(t N )dZ 3 ..-dZ N , (2.8) 

which is called a multilinear Fourier multiplier with symbol B, thus we identify the 
symbol and the operator. Whenever we write a symbol in the variables (£, £i, . . . £/v), 
it should be understood as a funciton of (£i, . . . £at) by substitution £ = + - • - + 6v- 
Hence for example, a product of single multipliers can be written as 

B o (0Si(ei)S2(6)[/,0] = BbHV){(Bi(-iV)/) • (fl 2 (-iV)</)}. (2.9) 

In addition, we introduce the following convention for bilinear multipliers. For 
any bilinear symbol i?(£ 1 ,£ 2 ) ) we assign the variables i] and ( such that i] = £ x and 
C = £2, but regarding (£,77) and (£, C) respectively as the independent variables. 
Hence the partial derivatives of the symbol B in each coordinates are given by 

(vfB, V V B) = (V 6 flfa, t ~ v), (V €l - V 6 )flfa, £ - 77)), 

(v^b, v c s) = (v €l s(e - co, (v 6 - v 6 )s(e - c, o). 

For any number or vector a, we denote 

(a) := J2+ lal 2 , a := C/(a) := ]4, #(a) := H(a), (2.11) 

|a| (a) 

which will mostly appear in the Fourier spaces, in particular U = a/— A/ (2 — A) 
and if = a/— A(2 — A). For any complex-valued function /, we often denote the 
complex conjugate by 

/ + :=/, /":=/, (2-12) 
to treat them in a symmetric way. Also we denote 

Jf — e~ itH xe itH f, f = Fe itH f, f ± = Te ±ltH f ± . (2.13) 

Remark that these operations are time dependent. 

Finally, we collect a few basic estimates, which will be used throughout the paper. 
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Lemma 2.1. Let 2 < p < oo, < 9 < 1, s G and a = 1/2 — 1/p. Then we have 



He-^||^ 2 < |tr ( ^||^- 2+3 ^(V) 2e V||^ 2 , (2.14) 
where p' = p/(p — 1) is i/ie Holder conjugate. For 2 < p < oo, we have also 



-^^|Up, 2 < |t|-^||^- 2+3 ^(V) 2e V|| L y, 2 . (2.15) 



The case 9 = is similar to the Schrodinger equation, but we gain U^' 2 ^ around 
£ = 0. The case 6* = 1 with |£| < 1 is the same as the wave equation, where we have 
worse decay for t — > oo, but better decay for £ — > 0, than the case # = 0. This reflects 
the fact that the group velocity for e~ ltH is bounded from below |V-ff"(£)| > y/2. For 
|£| > 1, the estimate with 9 > follows from interpolation between # = and the 
Sobolev embedding, hence it is just the same as for Schrodinger. 

Proof. The estimate in the Besov space is proved by the standard stationary phase 
argument after the Littlewood-Paley decomposition, and its interpolation with L 2 
conservation. See [T7] for the details in the case = 0. Then we can replace Besov 
with LP spaces by the embedding 

S° 2 CI7, B° p/j2 DL/, (2.16) 

for 2 < p < oo. The Lorentz version follows from real interpolation. □ 

From the above decay estimate, we get the Strichartz estimate (cf. [IT]). 

\\ e ~ lt fWL^mnu^L^B^ ^ 11/11 a-V6iABi 2+i i#i, (2-17) 

where we gain U 1 ^ for the endpoint norm. 

For bilinear Fourier multipliers, we will use the Coifman- Meyer bilinear multiplier 
estimate [7] (see also [H]): if 

1^5(6,6)1 < (161 + Nr H - |/31 , (2.18) 

for up to sufficiently large a, (3, then for any pj 6 (1, oo) satisfying l/po = l/p\ + l/p2, 
we have 

H-B[/> <7]||£P0(R<i) < ||/||LPl(Rd)||5f|| iP 2(Rd). (2-19) 

As was shown in [TJ], one cannot generally replace the right hand side of (12.181) by 
I6I j but we can reduce our regular multipliers to the above case. However 
the above estimate is not applicable when the multipliers are really singular due to 
divisors coming from the phase integrations in the main estimates for Theorem 11.11 
In this case, we will use another bilinear estimate (Lemma llO.ll) . which allows general 
multipliers with singularities. As for regular multipliers, the following specific form 
will be often sufficient: 

Lemma 2.2. Let d e N and k G N. Then we have 

I /. \ 2fc(l— a) / (■ \2ka 



sup 

0<a<l 



((6,6)) 



2k 



if, 9} 



^ ll/llx£i(R«i)llff \\l p x 2 (ri)i (2.20) 
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for any Po,Pi,P2 £ (1,°°) satisfying 1/po = 1/pi + I/P2, where (^1,^2) denotes the 
2d dimensional vector and so 

((£i,6)> = v / 2 + I6I 2 + I6I 2 - (2.21) 

Proof. Denote the above operators by B a , which is analytic on the strip < Re a < 1 
as bilinear operators (H d ) 2 — > L 2 . When a — or a — 1, B a satisfies (I2.18p . hence 
the estimate follows from Coifman-Meyer. Since the Fourier multiplier 
is bounded on IP uniformly for t G K, we get 

sup (t) 1+d (\\B it [f,g)\\ LP0 + \\B 1+it [f,g]\\ 

ten v ' 

which is extended to the whole strip < Re a < 1 by the three lines theorem for 
analytic functions. □ 



3. Normal forms 

The idea of using normal forms for asymptotic behavior of nonlinear dispersive 
equations goes back to |30j in the case of nonlinear Klein-Gordon equations. For the 
NLS, quadratic terms have in general resonance sets with one co dimension (except 
for u 2 ), which give rise to singularities in the normal forms. Here the resonance 
set is given in the spatial frequency where the time oscillation of the nonlinearity 
coincides with the linear evolution. We call it temporal resonance. 

In [Tj5] for the scattering of NLS in three dimensions, such singularity was avoided 
by applying the normal form to remove only some of the quadratic terms, leaving the 
others with better decay properties. Here the decay for the remainders comes from 
the fact that the oscillation phase of the nonlinearity relative to the linear evolution 
is not stationary at those spatial frequencies. We call it spatial resonance. When we 
simply say "non-resonant", it means that it is not simultaneously temporally and 
spatially resonant. 

Our transforms follow the same idea as [19J. However, the actual forms and 
computations are totally different from the NLS, due to the lower order term 2 Rev, 
and the special structure of the nonlinearity. In addition, we can not completely 
avoid singularities because of some degeneracy in the non-resonance property, for 
which we need another idea (a bilinear estimate for irregular multipliers). 

In the previous paper [18], we introduced the normal form M(u) as defined in 
(I1.15p . As is indicated in Theorem 11.31 this transform is quite natural in view of 
the nonlinear energy, and will be essential for the proof of Theorem 11.21 However, it 
does not fit the weighted estimate for Theorem 11.11 very well, where it suffers from 
a derivative loss. To avoid assuming higher regularity, we introduce another normal 
form which is more suited to the weighted estimate. For that purpose, we start with 
general bilinear transforms. 

First we rewrite the equation for u = ip — 1 = u\+iu 2 in terms of the real variables: 

ui = —Au 2 + 2uiu 2 + \u\ 2 u 2 , 

u 2 = —(2 — A)ui — ?>u\ — u\ — \u^U\. 



(3.3) 
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We consider a nonlinear transform given in the forrrJl 

w := u + B[[ui,ui] + B' 2 [u 2 ,u 2 ] (3.2) 

where Bj are real- valued symmetric bilinear Fourier multipliers, which we will specify 
later. From the above equation for u, we derive 

(id t + A-2Re)B' 1 [u 1 ,ui] 

= -(2- A)B[[u 1 ,u 1 ] + 2iB[[u 1 ,-Au 2 + 2u 1 u 2 + \u\ 2 u 2 }, 
(id t + A-2Re)B' 2 [u 2 ,u 2 ] 

= -(2 - A)B' 2 [u 2 ,u 2 ] +2iB' 2 [u 2 , (A - 2)m - 3m? - u\ - |w| 2 wi], 
and hence we get the equation for w. 

(id t + A - 2 Re)w = B' 3 [u 1 , ui] + B' 4 [u 2 , u 2 ] + iB' 5 [ui, u 2 \ 

+ |w| 2 wi +iC' 3 [ui,ui,u 2 ] + iC A [u 2 ,u 2 ,u 2 ] +iQi(u), 
where Bj (j = 3, 4, 5) are bilinear multipliers defined by, 

B' 3 = 3-(0 2 B[, B', = 1-(0 2 B' 2 , B' 5 = 2 + 2\&\ 2 B[-2(Z 1 ) 2 B' 2 , (3.5) 
C'j (j = 3, 4) are cubic multipliers defined by 

C 3 (Zi, 6, 6) = i + 451 (6, 6 + 6) - 6£? 2 (6 + 6, 6), 
Cl(£i,£ 2 ,£ 3 ) = 1-2^(6 + 6,6), 



(3.4) 



(3.6) 



and Qi(w) = 2_B^[mx, |m| 2 U2] — 2£> 2 [ii 2 , |^| 2, "i] gathers the quartic terms. 

Now we choose our B[ and B' 2 . Recall the diagonalized equation for v = Ui + iUu 2 : 

iv t — Hv = M v {u) := U (3m 2 + w 2 , + I m | 2m i) + i{2u\U 2 + |w| 2 -u 2 ), (3.7) 

where we observe that we have U~ l in each u 2 compared with v , and also on the 
imaginary part relative to the real part. To counteract this effect, we want 

£i = 0(|6li6i), B' 5 = 0(m 2 \), = 0(161 + 161 + 161)- (3.8) 

It is easy to see that B' A = B' 5 = is impossible due to the symmetry of B[ and B' 2 . 
The simplest choice for B' 4 = is given by 

B[ = B' 2 = (0~\ B' 4 = 0, B' 5 = A(£)- 2 ti 2 , C' 3 = C' 4 = U 2 , (3.9) 

where ££2 = £'6 denotes the inner product in M 3 . This leads to M(u) in (11.15p . and 
is the most natural from the energy viewpoint (11.211) . Moreover, the cubic part has 
a subtle but remarkable non-resonance property, which will be revealed in Section 
H] for the proof of Theorem 11.21 

On the other hand, the simplest choice for B' 5 = is given by 



-^ = ^ = ((6, 6)) 



-2 1 



2 + I6P + I6I 



4 "2 + |6| 2 + |6| 2 ' ^ 5 ~°' 64 



16 + 61 


2 


+ 161 


2 


2 + 16 + 6! 


2 + l6l 2 



(3.10) 



We will introduce Bj and Cj afterward by transforming £?'• and C' . 
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where we omit C 3 since it is just complicated without any important structure. This 
choice is not directly linked to the energy, and the cancellation in the cubic part 
is weaker. Its advantage is that we do not lose a derivative for (£) <C |6| ~ |£ 2 | 
as in B' 5 of (13. 9p . Such a derivative loss for high-high interactions can be saved 
by increasing the regularity. However, in the weighted estimate for bilinear terms, 
we will encounter another derivative loss (I10.33[) by integration on the phase, and 
therefore we would need H 2 to close the estimates with the choice (13.91) . 

The latter derivative loss can be observed in the Schrodinger case as well: look at 
the weighted estimate for the typical bilinear term |e~ l * A y?| 2 in the Fourier space 

^ e «A| e -*/^| a = J e «n^(_^)^ - v )d V , (3.11) 
where Q := |£| 2 — M 2 + |£ — v\ 2 - Partial integration on e ltQ gives a factor of the form 

(3.12) 



|V^i_|2£-^| 



|v„fi| iei 

which causes a derivative loss in the region |£| <C \rj\ ~ \rj — £|- However, those 
two normal forms eventually lead to the same asymptotic behavior, in the case of 
weighted energy scattering. See Section ISTTl 

We remark that there are other possible choices. For example, if we want decay 
in B' 3 as well, we may choose 



B[ = (0~ 2 (3 - 4<£)" 3 £i&), B' 2 = (0~ 2 (l - 4<£>" 2 £i6) 

>-2, 



B' 3 = B' A = 4(0"^i6, (3-13) 

B' 5 = 4(0~ 4 (4^o6 + 4|e | 2 |6| 2 + 2(£i6) [£„(&> + 46 - 4f 2 )], 

but it suffers from derivative loss in all B'j, so we will not use it. In the rest of the 

paper, those notations B'-, C'- and Qi are reserved for the second choice (13.101) . 

Now we fix the notation for our transforms corresponding to (13.91) and (I3.10p 
respectively b}0 

z : = m{u) = v + (er 2 N 2 , , . 

-2 -2 ( 3 - 14 ) 

Z:=v + b(u) :=v- ((6,6)) «i] + ((6, 6)) ["2, "2]- 
Then we get coupled equations for (t> , 2) and (t> , Z) in the form 

v = z ~ (£)~ 2 \u\ 2 , iz-Hz = N z {u), (3.15) 

v = Z-b(u), iZ - HZ =N z (v). (3.16) 
The nonlinear terms are given by 

M z {u) := U {2u\ + MV} - i{V)~ 2 V ■ {4u!Vu 2 + V(\u\ 2 u 2 )} , (3.17) 

Mz{v) := B^v^yx] + B 4 [v 2 ,v 2 ] + C^v^v^Vi] + C 2 [v 2 ,v 2 ,v 1 ] 
+ iC 3 [v 1 ,v 1 ,v 2 ] + iC 4 [v 2 ,v 2 ,v 2 ) + iQi{u), 



(3.18) 



ir The pair (v, z) in the previous paper |18j is equal to U 1 (v, z) in this paper. 
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with the bilinear multipliers Bj defined by using Bj in the case ( I3.10p . 




(3.19) 



and the cubic multipliers Cj defined by using Cj in ( 13. 6ft . 



1 




(3.20) 



Q 1 H = -2((ei,6)) _2 K 



u\ 2 u 2 ] - 2((£i,£ 2 )) 2 K, M 2 wi]. 



(3.21) 




(3.22) 



i 



More precisely, they are written as linear combinations of products of those terms 
on the r.h.s. of (13.221) and regular multipliers as in (I2.20p . together with the Riesz 
operator. 

Remark 3.1. By the transform u i— > v := u — iH' 1 ^, the Boussinesq equation (II. 3p 
can be rewritten as 



hence it can be regarded as a simplified version of our equation for (v,z). In fact, 
our proof for Theorem 11.11 directly applies to this equation, while Theorem 11.21 and 
Theorem 11.31 rely much more on the nonlinear energy structure, which is different 
for (II. 3p . For the (physically relevant) two dimensional case, refer to [18] for a final 
data result, which is also applicable to this equation. 



In this section, we consider the final state problem in the energy space, proving 
Theorems ll.2l and ll.31 Our strategy for Theorem ll.2l is essentially the same as [25] for 
NLS, using compactness to get weak convergence, and energy conservation to make 
it strong. The normal form transform M{u) is essential both for the uniform bound 
and for the equicontinuity of approximating sequence of solutions. In addition, we 
will use a hidden non-resonance property in the cubic terms of ( 13.151) to get the 
equicontinuity. It is noteworthy that the argument in [28] does not work for NLS 
without L 2 finiteness, whereas for GP we are using only the energy conservation. 
We start with the proof of Theorem 11.31 since the result will be used in the proof of 
Theorem 11.21 

Proof of Theorem \1.3[ The identity (ll.2ip is straightforward. We construct R by 
solving the following equation for g with any given / G B X {L } k): 



iv t + Hv = U{v{), 



(3.23) 



4. Final data problem 



9l = R f ( 9l ) := h - (V)- 2 {g\ + {U- l f 2 ) 2 } , g 2 = IT 1 /* 



(4.1) 
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We use the following product estimate for general u,v: 

ll<Vr 2 H|£x< \\U{uv)\\v ~ \Muv)\\ B -x 

< WvWhiWVuWh-^ + IMU°°+z6||Vt>||i,2 

< 



(4.2) 



(4.3) 



(4.4) 



\\Uu\\» + \\u\\)/ e 2 \\Uu\\)/ 2 

Applying it to Rf(gi), we get 

IKV)- 2 ^ 2 2 || H1 < ||^ 2 2 || L2 < ||/||^||/||^ 2 , 

\\ R f( u ) -Rf(v)\\ H i < Wu + vW^Wu + vW^eWu-vWfji, 

and in addition, 

\\^r 2 9 2 Am < \\Ufi\\» < lkiV^|| L3/2 < IbilUall^H^, 

ll(V)- 2 ^|Ua< ||s 2 || L 3 = Mi 6 . 

Hence the map g i— > Rf(gi) + ig 2 is a contraction on the set 

{g E F 1 {2L,C{k + L^k 1 I 2 )) \ g 2 = U^fr}, (4.5) 

for the H 1 norm and some C > 0, if k > is sufficiently small. Hence it has a 
unique fixed point g\ = Rf(gi), by which we can define R(f) = gi + iU~ 1 f 2 - Then 
we have M o R = id and the uniqueness of R follows from that for the fixed point. 
For the continuity of R, let /, g G H 1 ^, k) and u = R{f), v = R(g). Then we have 
by (53), 

\\U(\u\ 2 - M 2 )|| L2 < Wu + vWfju + vW^Wu-vW^ « \\u-v\\ H i, (4.6) 
if K is small. Moreover, 

\\ u - v \\ 2 m = IK-^ill^i + \\h-92\\m- (4-7) 
For the real component, we have 

hi-vi\\Hi = \\RfM-Rg(v 1 )\\ A i 

< ll/i-^ill^ + IIW^K-^ll^ + IKV)- 2 ^-^) 2 -^- 1 ^) 2 !!^- 

where the second term is bounded by 

IK + ui|Ue|K - v i\\m < IK - viWiPi (4-9) 
and the third term is bounded by using (14.21) . 

||/ 2 -^||^||/ 2 + ^t / xll/ 2 + ^|li / 6 2 « \\f2-g2Wm. (4-10) 
Hence for any e > 0, we can get by choosing k > small enough, 

IK - vxWhi < \\fi - gi \\ m + e\\f 2 - g 2 \\m, (4.11) 

and the desired Lipschitz condition follows. □ 
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Proof of Theorem \l.S\ First we rewrite the equation ( 13. 15ft for z = M(u), using the 
renormalized charge density q(u) := \tp\ 2 — 1 = 2u\ + \u\ 2 \ 

iz t - Hz = Af z \u)+Af 2 (u), 

Ml := ^(q 2 - \u\ 4 - 2u\) - • {2<?Vw 2 + m 1 (2m 2 Vm 1 - «iVw 2 )} , (4.12) 
1 (V) 

A^ 2 := -U{u\u x ) - i(Vy 2 V ■ {u\Vu 2 ). 

Let L = 2||(V)2; + ||l2 + 1 and z° := e~ ltH z + . Since ||2; (t)||i6 — > as t — > oo, for large 
T > we have 2°(T) e H X {L, k). Let V ;T be the global solution of (jl.ip satisfying 

V^ T = 1 + « T , u T (T) = R{z°{T)), (4.13) 

and let z T = M(u T ). Then we have z T (T) = z°(T) and the above equation (I4.12p 
with (z,u) = (z T ,u T ) and q = q T = q(u T ). Moreover we have 

||(V)^|| L2 < E^fl 2 = S(R(z°(T)),0) < 2\\(V)z + \\ L2 , (4.14) 

so z T , u T and q T are uniformly bounded in H 1 , Fx and L 2 respectively. By the 
Sobolev embedding, u T is bounded in L\. Combined with the L 2 bound on q T and 
the fact that (L 3 + L 6 ) n L 2 C L 3 , it implies also that and g T are bounded in L 3 . 

Moreover, u T and z T are equicontinuous in C(R; iS'(R 3 )) by their equations. Hence 
for some sequence T — ► oo, u T and z T converge in CiM.^w-H 1 ) and C^w-if 1 ) 
respectively. We denote their limits by u°° and z°°. Then g T also converges to q°° = 
q(u°°) in C(R; L 2 loc ) and C(R; w-L 2 ). Also it is easy to check that := + illu™ 
and z°° satisfy the equation (13. 15p . 

To get weak convergence of e ltH z T globally in time, we need its equicontinuity at 
t = oo, i.e., for any test function (p e iS(M 3 ) we have 

t2 (p\e UH Af z (u T (t))) x dt^O (4.15) 

L 

as £ 2 > t\ —* oo uniformly in T. Those terms in N}{u) are bounded respectively by 
(here we omit T), 

lk 2 ||z,12/ll ^ IMIi2 nL 3, 

II^|w| 4 |Il^/ii < \\u 3 Vu\\ L i < ||n|| 3 6 ||Vn|| L2 , 

\\U(ul)\\ L u/u < ||u?Vui|| L i < ||wi||is n L«l|Vu||ia. (4-16) 
||gVw 2 |Ui2/n < ||g||L2nL3||Vu|| L 2, 
||?iinVu|| L i2/n < ||ni|| L 3 nL 6||M|| L 6||Vn|| L 2, 

hence by using the L 12 decay of e ltH (" 12. 15j) . 

T \{ V \j tH N x z {u T {t))) x \dt < ft- 6 '** < t~ 1/4 -> 0. (4.17) 

The remaining term A/" 2 is bounded only in Lj , which gives by the L p decay the 
critical decay rate 1/t, and so we cannot get compactness in this way. More precisely, 
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when we decompose them into dyadic frequencies by (12.61) 

A/? = - J2 ^{x(V) i (^)x(V)S} + i^ ■ WV)1«l)x(V)%} , (4.18) 

we have trouble only when j <C min(l, k), since otherwise U or V outside gives the 
factor j and so those terms are bounded in L 12 ^ 11 as above. In order to bound the 
remaining part, we use a non-resonance property due to its special structure: 

Ml = ~U(u 2 2 v) - i(V)- 2 {H(£)U(b) - £6} [ul u 2 ], (4.19) 

where the second bilinear multiplier has the small factor £ 1; since it equals 

(0" 2 {iei((e)^(6)-iei) + ier-a 2 } 

leil&l Ki l (a +26) -6 , ( 4 -2o) 



+ 



«0 + <&»<&> + J (0 2 <0 2 ' 

and the multiplier inside the big braces is bounded by the Coifman-Meyer estimate 
(I2.19P if restricted to <C |£ 2 | ~ So this term is bounded in L^ 11 . 
For the first term in (14.191) . we integrate on e tt<yH ^ +H ^ 2 ^ 

" h e u ^ +H ^U[u 2 2 ,^}dt 



[-ie itH 'B[u 2 2 ,v]]% + [ 2 ie itH B[2u 2 ii 2 ,v] - e itH B[ulK(u)]dt, 



(4.21) 



with the bilinear multiplier B defined by 

b = + H^-'uio = (o- 2 {l + Hi&H^rT 1 . ( 4 - 22 ) 

where the multiplier in the braces is bounded by Coifman-Meyer (I2.19P if <C 
|£ 2 | ~ From the equation (13. ip . we have 

u 2 = -q + A Ul - q Ul E L?H-\ 

M(u)=U(\u\ 2 + qu 1 )+tqu 2 eL™(L 3 / 2 nL 2 ) x , 

so u 2 U2 e ij 1 x /f^ 1 is bounded in H^ 1,3/2 fl #~ 2 . Hence ( 14.21ft restricted to 
|6| < min(l, |6|) is bounded in t^ 7 ^ + h 1/4 L 12 . 

Thus we conclude that z T := e ltH z T is equicontinuous in C([— oo, oo]; iS>'(IR 3 )), and 
so by choosing appropriate sequence T — > oo, we get 

i T -> i°° in C([-oo, oo]; w-tf 1 ), (4.24) 

together with i°°(oo) = z + from the initial condition z T (T) = e~ tTH z + . By the 
lower semi- continuity for T — > oo (along the sequence), we have 

Ei(i>°°) = WW^Wh + h\U\u°°\ 2 \\h < liminf E^ T (T)) = \\(V)z + \\ 2 L >, (4.25) 

Z T—>oo 

and for t — > oo, 

\\(V)z + \\ 2 L2 <limmfE 1 (i>°°(t)) = E 1 (i/;°°). ( 4 .26) 

t— >oo v ' 

Hence both inequalities must be equality, which implies that the convergence is 
strong in both cases. Hence e ttH z°°(t) — > z + strongly in H\ as t — > oo. □ 
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Remark 4.1. We can show that for some asymptotic profiles z + G H 1 , the above 
constructed solution has V\ unbounded in L 2 X for large t (v 2 is bounded in L 2 by the 
energy). So we really need the nonlinear map M to have the scattering statement 
( I1.16P in H 1 . The proof goes by contradiction. If v% = U\ is uniformly bounded in 
L 2 , by the L 2 bound of q, \u\ 2 is also bounded. Then it is easy to see that M z {u) is 
bounded in UL P X for 1 < p < 6/5. Thus we have by the LP decay with 1/p := 1/6 — s 
for e > small, 



T 

< f S^\\KHs))\\ ULl/ e +s ds <T~ 3£ , 



and so by letting t — > oo, we get 



z+ = te tH/k k 3 / 2+3 "f(kx), 

fcS4- N 



(4.27) 



e~ iTH z + - z(T) e U 4/3+£ L p x . (4.28) 



Since lJ- l z 2 {T) = u 2 (T) G L^L X C H-^ 2 ~ £ / 3 ' p for T large, we have 

(Ime- 4TH z+) <1 = (C/« 2 (T) + ?7 4 / 3+e L5?Lg)<i C U A/3 L^L P X , (4.29) 
which is clearly not implied by the assumption z + G iJ 1 . For example, let 



(4.30) 



for some / G S satisfying supp / c{l<|£|<2}. Then z + G H 1 but 

WU-^-iH/k^ > ^-4/3+3/2+3^-3^1^11^ ^ ^ (i/^^). (4 . 31) 



However it is not clear whether ||fi(t)||z, 2 is growing, or infinite for all time (v (t) G 
L 2 X persists in time by the equation for v). We expect that both cases do happen, 
depending on the asymptotic data. 



5. Initial data problem 

The strategy for Theorem 11.11 largely follows that in [19] for the NLS: we will 
derive an a priori bound on the weighted energy norm of the solution after removing 
the free propagator, namely e ltH v{t). Then the a priori bound gives space and time 
decay properties of the solution by getting back the propagator, and this decay 
property is used in turn to bound the Duhamel integral for the a priori estimate. 
However, the linear decay property is not sufficient to bound the quadratic terms, so 
that we have to integrate by parts on the phase to gain more integrability and decay 
from the oscillations. At this point we meet the special difficulty of our equation, 
that is the singularity due to degeneration of the nonresonance property around 
£ — > 0. To treat it, we will derive a bilinear estimate for singular Fourier multipliers, 
and investigate the shape of the degeneration in a similar way as in our previous 
paper [TS] for the two dimensional case. 
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Now we set up our function spaces. We will derive iterative estimates for both v 
and Z of (13.161) in the following space-time norms. Recall the notation (12.1 3p for J. 

\\Z(t)\\ x{t) := \\Z(t)\\ m + || JZ(t)\\ H i, \\Z\\ X := sup \\Z(t)\\ x{t) , 

* (5.1) 

\\Z\\s ■= \\Z\\ L f>Hl + \\U 1/6 Z\\ L 2 H i,6. 

The last norm is finite for e~ itH ip if <p E H^M. 3 ), by the Strichartz estimate f[27LTl) . 
The key ingredient is the H l bound on JZ, which can be written in the Fourier 

space 

\\JZ\\ Hl ~ IKOV^I^. (5.2) 

More precisely, we are going to prove 

ft 



e 

T 



i{t - s)H M z {s)d,i 



< (T)- £ , (5.3) 

X(T,oo) 

for some small e > 0. This estimate implies the desired scattering for Z in (x)~ H 1 . 
The scattering for v is the same, because we will prove that the difference b(u) is 
vanishing faster in the same space X. 

5.1. Decay property by the weighted norm. We derive a few decay properties 
as t — > oo and £ — ► of the X space (15. ip . First, the commutator relations [Vj, Jk] = 
5 jtk and [(V), J] = -(VrV imply 

\\JVv(t)\\ L 2 + \\J{V)v(t)\\ L , < \\Jv(t)\\ H i + \\v(t)\\ Ll . (5.4) 
Since l/\x\ G L '°°(R ), we have by the Sobolev and the Holder, 

< ||e^(t)|| L s /5 , 2 < \\xe itH v(t)\\Ll ~ \\Jv(t)h>. { ' j 

Thus we obtain 

\\U~ 2 v\\l% < \\v\\h-i + \Hm x ~ H^Hfli < \\v(t)\\ x{ t). (5.6) 
Choosing 6 = and p = 6 in the L p decay estimate (I2.15p . we have 
\\(V)U-^v(t)U, < t- l \\(y)e itH v(t)\\ L l^ 

Combining with the Sobolev bound (15. 6p . we get 

|||V|- 2+5e /Vi(t)IUa < mm(l } t- 9 )\\v(t)\\ x[th 

IllviViWIUg < minir^r^Ht)^, 

for < 9 < 1, where t><i and t>>i denote the smooth separation of frequency (12.71) . 
In particular, we have 

I|tf-Mt)llz5 £ <*>~ 3/5 IM*)ll*M. ll^llx?^ £ IMIsnx- (5.9) 
Also we obtain the Strichartz bound on u = V\ + iU~ 1 v 2 : 

\\u\\l°°w ^ H^'MIl-h 1 < IMU, 

, (5.10) 

\\u\\l^H^ < F>1 IW' 8 + \\U W<i|| L 2L6 < v xns- 



SCATTERING FOR GROSS-PITAEVSKII IN 3D 19 

For the estimate on b(u), it is better to use the wave-type decay. Choosing p = 4 
and 9 = 2/3 in (|2.15p . and using complex interpolation, we have 

||<v> 2/3 fryt)|| L 4 < t- 7 i l2 \\u- l /\vy tH v{t)\\ T ^ 



L 

^r^lKvje^wil^llcr^vje^wil^ (5.11) 
<t- 7 '"\\v(t)\\x(ty 



6. Estimates on the normal form 

In this section, we derive decay estimates on b(u) (defined in (13.141) ). and the 
invertibility of the mapping v i— > Z. First from the LP decay property of v and the 
bilinear estimate (12.201) . we have 

\\b(u)\\ H >, P < licr^n^iitr^llip < licr^ii^ncr^llig 

< <r 3 * /5 IK*)lll( t ), 

for < < 2, where 



10 11 . , 

- = 1-77 = - + -, 2< Pl ,p 2 <6. 6.2 
p 3 pi p 2 

For the weighted norm, we compute Jb(u) in the Fourier space, using the notation 
in (I2.13p . It is given by a linear combination of terms of the form 



e -ftfftt) V{ /" e^ )TH( " )T//(c) fi(r/,C)^ ± (r/)t; ± (C)^ 

= e -itH(0 j e itn ^y(v) B + , tV f } fi • 5)^(77)^(0 + ^ ± (^)Vt) ± (C)]^ ( 6 - 3 ) 

= .f UvfB + itvfn ■ b)^, u ± ] + (j w ) ± ]j , 

with B = ((6,6)>~ a or ((^k))- 2 ^!)- 1 ^)- 1 , and 



n:=H(0TH(r))TH((), V^O = Vif(£) T Vif(C). (6-4) 

By the bilinear estimate (I2.20p and the L p decay property (15. 6p . (15.91) and (15.111) . 

we have 

||Vf B[v±,v±]\\ Hi < \\U^v\\ L3 J\U- 2 v\\ Ll < (ty 3/10 \\vr x , 

\\tvfn-B[v±,v±]\\ Hi < \\t^u- l v\\l % < (tr 1/6 \H x , (6-5) 
(JtOlIk < wu-MlLiWu-'Jvh, < (t)- 3/10 \\v\\ x . 

Thus we obtain 

\\b(u)(t)\\ Hi < (tr 9/1 °\\v\\ 2 x{t) , \\Jb(u)(t)\\ m < (C^IMIW 

\\Ku)\\s(T u T2) £ (Tir 7/1 °\\v\\ x{TuT2) , 

for any t > and any T 2 > T\ > 0. Hence by the Banach fixed point theorem, the 
map i)hd + b(u) is bi-Lipschitz in small balls around of X(t) uniformly for all 
t el, and also globally in X D S. 



(6.6) 
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6.1. Other normal forms. It is a natural question if we get the same asymptotics 
for other normal forms, in particular, for z defined in (13.141) . It is easy to see that 
the above argument works to get 

n<vr 2 M 2 iu w < <*r 1/6 iM&. (6.7) 

For this, we need to modify the second inequality in (I6.5P for the interaction with 
|£| + 1 <C |£i| ~ | £2 1 because we lose regularity. Using the Schrodinger type decay 
estimate in that region, we get nevertheless 

HfV^n • (erV^llk £ P 1/2f/ ~Mlit + ll^>ill^lk>ilUi . 8) 

— 1/6 11 i|2 \ ■ ) 

% (t) \Mx- 

Therefore we have 

Mt)-v(t)\\x (t) < (t)- 1/6 \\vr x , (6.9) 

in other words, z, v and Z have the same asymptotics in the strong topology of X, 
as long as v is in X. 

7. Initial condition 



Here we check that the smallness condition f 1 1 . 1 2 j) on u(0) is equivalent to smallness 
of v(0) inX(O). 

By Sobolev and Holder in Lorentz spaces, (11.121) implies 

IK0)IUi < ||V W (0)|| 4 / S , 2 < \\xVu(0)\\ Ll < S. (7.1) 

Hence we have 

||u(0)||fli + \\x Ul (0)\\ Ll + \\xV Ul (0)\\ L 2 + \\xVu 2 (0)\\l* < 5. (7.2) 
Moreover, by using the commutator [x, U] = JF _1 [iV^, f/(£)]JF = iV^U, 
\\{V)xUu 2 \\ Ll < \\(V)Uxu 2 \\ L 2 + |KV)(V e C0u2||Li 



< ||Vxu 2 ||li + \\u2hi: 



(7.3) 



hence we get || (x)v (0)||#i < 8. 
On the other hand, since 



[Jj, V k U- 1 } = -5 hk U- 1 + V^iVfcET 1 , (7.4) 

we have 

\\JVu\\ L% + ||VJu|U| < \\U- l VJv\\ L i + \\[J,VU- l ]v\\ Ll + \\u\\ L i 

< \\Jv\\ Hl + \\u- l v\\ L% < IKt)!U( t) . 

In particular, by letting t = 0, 

\\(x) Ul (0)\\ L * + [|<ar>VT*(0)||xa S IK*>«(0)lk- (7-6) 
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8. Estimates without weight 

In this section, we estimate Z in the Strichartz norm LfB.\ H U 1 l % L 2 t Hlf' ^ by 
putting the nonlinear terms Afz(v) in the dual Strichartz space, as in (12.171) . 
The bilinear terms are bounded by using (12.201) : 

\\ B A v± i v± ]\\ L 4 t / 3 Hi> 3/2 ~ NU°°.ffi|MUvsii6 < \Mxns> (8-i) 

for j = 3, 4, where the last factor is bounded by the decay estimate (15. 8p . On interval 
t G (T, oo) for T > 1, we can estimate instead 



|2 



< 1/2|| 7 .|| 1/2 |U 7 .||3/2 < T -l/2|. , 

~ J II ^IIloo^I ||t-f || 1,6 ^ J H^llX- 



(8.2) 



The cubic terms Cj with j = 1, 2, 3, 4 are similarly bounded by 

\\C j [v ± ,v ± ,v ± }\\ L2tH i, 6 / 5 < \\U- x v\\ l ~hA\ U ~ 1v \\1±ls ^ Nix, (8.3) 



by using (15. 6ft and (15.81) . together with the regular bilinear estimate ( 12.201) . For 
t > T we gain at least 7" 2 (- 3 / 5 + 1 / 4 ) = y-7/io f rom ^ e jj^ decay of U~ 1 v. The quartic 
term Qi is estimated by using (I5.10p . 

II Ql(«) |L*/8ffi,3/a < ||u||i4Hi,3||u|U<»L6 < |MlW (8.4) 

For t > T we gain at least T~ 3 / 5 from the decay of it = v± + iU~ 1 v 2 . 
In conclusion, we have 



e -i(t-s)H N ■ dg 



T 



^ Wz\\lZ >T Hi,*/*+Ll >T Hi 
S(T,oo) (8.5) 

<^ 1/2 (IMIxn5+ IMIlns), 
so once we get a uniform bound on \\v \\xns, the scattering of Z and t> in H 1 follows: 

3!^ + G if 1 , \\Z(t) - e- itH v + \\ m + \\v(t) - e- itH v + \\ m < (t)~ 1/2 , (8.6) 

where we used ( 16. 6 j) . 

9. Estimates without phase derivative 

Now we proceed to the LfH]. bound on JZ. First we rewrite the equation for Z 
by replacing v with Z in the bilinear terms: 

5 

Mz{v) = B 3 [Z U Z x ] + B 4 [Z 2 , Z 2 ] +J2 C i + Q^ u ) + Q^ u ) =■ Z ^ ( \ 

C 5 (v,v,Z) := -2B 3 [b(u),Z 1 ], Q 2 (u) := B 3 [b(u),b(u)\, 

because we are going to integrate on the phase in time (in some Fourier region), 
where the difference of oscillation between Z and v will become essential. In the 
higher order terms, it is negligible thanks to the time decay of b(u). 
Applying J to Af' z , we get bilinear terms in the Fourier space like 

jf J Vf {e^Bjirj, £ - ^(s, V )Z±(s, £ - v ) } d V ds, (9.2) 
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with j — 3,4 and the phase Q := H(£) =F H(r)) =F H(£ — rf). If the derivative 

(ri) 

lands on Bj, then since V) Bj is a bounded multiplier, its contribution is estimated 

t 



as above. If vl lands on Z*, we get terms of the form in the physical space 



f e^B^^JZ^ds, (9.3) 
Jo 

whose contribution in H\ is bounded by using the Strichartz together with the 
regular bilinear estimate (I2.20p . 

\\Bj[Z ± 1 (JZ) ± }\\ L iri H i,z/2 < 1 1 ^|| ^4/3^1,6 1| JZ || L oo H \ < \\Z\\ 2 Xf]S , (9.4) 

where the L 4 / 3 if 1,6 norm was bounded by the Strichartz 5* for small t and by the 
decay (15.81) for large t. Hence if we restrict on i 6 (T, oo), then we get additional 
factor T -1 / 4 . 

Similarly, we have cubic terms like 

V e T // e^C J (6,6,e3)^(s,ei)^(s,6)^ ± (5,e3)«2^, (9.5) 

for 1 < j < 4, where = #(0t#(6i)T#(6)T#(6)- If V ? hits its contribution 
is bounded in the same way as before. If Vg hits u , then its contribution is estimated 
by using the Strichartz and the bilinear estimate (I2.20p : 

\\C j [v ± ,V ± , (Jf) ± ]|| i 2 H i, 6 /5 < || Jv\\ L ^ H i\\U' l v\\ L 2 H ifi\\U' 1 v\\ IjCX>L & 

+ \\U~ l Jv\\ 
< \HxnS, 

where if the derivative in the H^' 6 ^ 5 norm lands on Jv (with large frequency), it is 
dominated by the first term on the right, otherwise we use the second term. Refer 
to (15.91) for the L 2 H l,& norm on U~ 1 v. For large t > T, we gain at least T~ 4 / 2 . 

For C5, we have just to replace the last v with Z in ( I9.5p . hence the final bound 
in (I9.6P is replaced by IMIinsll-^ll^nS- 

The quartic terms Qj(u) with j = 1,2 are regular. So we write in the physical 
space 

j f e -i(t-s)HQ.( u j ds = f e -i(t-s)H( x _ s ^H(C))Qj(u)ds. (9.7) 
Jo Jo 

Then by Strichartz, their contribution in L°°if 1 is bounded by 

1 1 XQj(u) 1^2^1,6/5 + \\tQj(u)\\ L 2 H 2, 6 /5. (9.8) 

For the first term we need to estimate xu. Since Ju = xu — t(V^H)u, we have by 
using (17. 5p . 

IMUg+tLg < \\M\in + \\(VzH)u\\ L 2 x < \\v{t)\\ x{t) . (9.9) 
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Then we have 

(9.10) 

and 

\\tQj(u)\\ L 2 H 2,6/5 < \\t 1/2 u\\ 2 LOOL(i \\u\\ L 2 L a\\u\\ LxL 3 < ||f||x nS . (9.11) 
For t > T, we gain at least T" 1 / 4 from the decay of u in L°°L 6 n L°°L 3 n L 2 L°° by 

d53D. 

Thus it remains only to estimate the terms with the derivative landing on the 
phase e lsQ in Bj and Cj. 

10. Estimates with phase derivative in bilinear terms 

This and the next sections make the heart of this paper. If the derivative lands 
on the phase in the bilinear terms, then the pointwise decay as above is at best 
t which is far from sufficient. So we should exploit the non-resonance property, 
through integration by parts in rj and s. We decompose Bj with j — 3,4 smoothly 
into two parts: 

Bj — Bf + Bj, (10.1) 

where Bj is supported in (£, rj) where the interaction is spatially non-resonant, and 
Bj in the temporally non-resonant region. Here non-resonance means simply that 
either 77 or s derivative of e does not vanish. 

We integrate the phase in 77 for Bj and in s for Bj. The strict intersection of 
the spatially resonant and temporally resonant regions is only at £ = 0, which can 
be compensated by the decay at £ = of the bilinear forms. This is the reason why 
the above decomposition is possible. 

The same type of argument has been used in [18J for the 2D final data problem. 
The main difference from there is that now we need L p -type bilinear estimates to 
close our argument for the initial data problem, and that we are not free to integrate 
by parts because the entry functions are the unknown solutions, whereas in the final 
data problem they were the given asymptotic profile. 

After the integration by parts, we get multipliers whose derivatives have much 
stronger singularity than allowed in the standard LP multiplier estimate, even if it 
were linear. The singularity is due to the behavior of H(£) around £ —>■ 0, which 
is close to the wave equation and thus enhancing resonance between the parallel 
interactions. Roughly speaking, the singularity increases for each derivative twice 
as fast as for the standard multipliers, even with the best decomposition of B x + B T . 

To overcome this difficulty, we take advantage of the room in decay in the 3D 
case to reduce differentiation of the symbol, employing a bilinear estimate with loss 
in the Holder exponent. It turns out that there is a narrow balance between the 
singularity and the loss in L p , such that we can close all the estimates. 

Specifically, we use only l/2 + £:or 1 + e derivatives for small frequencies, while 
the standard LP estimates need at least 3/2 + e derivatives of the symbol. Moreover, 
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it will be important for us to exploit the smallness of region where the singularity 
is the strongest. 

The plan of these two sections is as follows. First we prove the bilinear estimate 
allowing some singularities. Secondly we carry out the estimates on all the bilinear 
terms, assuming some bounds on the bilinear Fourier multipliers with divisors. Fi- 
nally in the next section, we prove those bounds for each multiplier, using geometric 
properties of the resonance sets. 

10.1. Singular bilinear multiplier with Strichartz. We introduce mixed (semi- 
) norms VB^ r for the symbols by 

\\f(^v)\\i/Bi r '■= lb'V(V)„/(£,77) Wqitit&PxR****), (10.2) 
where x J is as in (12.61) . 

Lemma 10.1. Let < s < d/2, and (p,q) be any dual Strichartz exponent except 
for the endpoint, namely 

l<p<2, Kq<2, - + - = 2 + -. (10.3) 

p q 2 

Let (pi,qi) and (p 2 ,<?2) satisfy 

111111 1 lis 

Pi P2 P qi q2 q g(s)' q(s) ' 2 d' (10.4) 
P<Pi,P2<oo, q<qi,q 2 <oa. 

Then for any bilinear Fourier multiplier B we have 
z itH B[u(t),v(t)]dt 



~ ^h^Bl^+L^B^J^L^LlAML^Ll^- (10.5) 



where the first norm of B is in the (£, rj) coordinates and the second in (£, () = 

(£,£-*?)■ 

Remark 10.2. q(s) is the Sobolev exponent for the embedding 5| x C L q<yS \ and 
l/q{s) gives the precise loss compared with the Holder inequality. 

The above estimate is a sort of composition of the Strichartz and the bilinear 
Fourier multiplier. Indeed, by choosing u(t,x) = f n (t)<p(x) and v(t,x) = ip(x) with 
11/nlU 1 < 1 and f n -^5e <S'(R), we can deduce 

Corollary 10.3. Let < s < d/2, and let (51,92) satisfy 

1111 ..lis , 

- + - = ^ + ^T> 2< qi ,q 2 <q(s), —-:=--- 10.6 
qi q 2 2 q(s) q{s) 2 d 

Then for any bilinear Fourier multiplier B we have 

\\B[<p,if>)\\i4 ~ \\ B \\L^B iltV+ L^B h Jv\\L%4^\\L^- (10.7) 
Practically, the above norm on B can be estimated by the interpolation 

ii^Hl 00 ^ ~ ii s iii«^. ii s iix«ff-i' ( io - 8 ) 

for s = (1 — 9)s + 9s\, with s 7^ s i an d 6 G (0, 1). For intermediate exponents, 
we can instead apply multilinear real interpolation to the above estimates, replacing 
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L°°i?2i by L 00 ^^ for free. But this does not work in the boundary cases, and 
indeed the estimates would become false, for example when s = or s = d/2, which 
will be used later. We abbreviate those norms by 

[H'\ := LfH* + LfH s [B s ] := LfB^ + LfBfa. (10.9) 

The proof given below essentially contains one for the corollary, but it seems 
difficult to decouple the lemma into the Strichartz and a time independent estimate. 
On the other hand, the corollary just barely fails to be sufficient for our application. 
In fact, more use of the Strichartz simply causes more loss in decay, compared with 
the weighted decay estimate. That is why we do not pursue the endpoint Strichartz, 
which corresponds to the worst case. Nevertheless, we prefer the flexibility in the 
space-time Holder exponents provided by Strichartz. 

Proof of Lemma \10.1\ By symmetry between rj and (, it suffices to prove for the 
norm in the (£, rj) coordinates. In the case s = 0, we get the above just by Plancherel. 
For s > 0, we want to integrate by parts in the Fourier space. Since the Fourier 
transform of u and v themselves are not smooth, we will use the cubic decomposition 
in x. The Fourier transform of the bilinear operator is given by 



10.10) 



X := J V )dri, /(£, V ) := J e itH ®u(t, V )v(t, £ - r,)dt 

= Hfy J e itHx [u(t, x + y)v{t, x)}dt. 
First we dyadically decompose B in the frequency for rj: 

X=^X n , X n := f B n fdr), B n := x n (V)„fl(7j,£ - V ). (i .n) 

nS2 z J 

By the definition of our norm on B, it suffices to bound each X n uniformly. Next, 
we decompose u(t, x) and v(t, x) into cubes of size n in x. Fix c e C^(R d ) such 
that for all i G K d we have 



^c(x-j) = 1. 



(10.12) 



Then for any ip G iS'(R d ) we have the following decomposition 



#0 = E ' ' ! '-'V;(0. V](x) := c(x/n - jMx), 

j&d (10.13) 

^(0 := e m ^](0 = J*[c(x/n)<p(x + nj)], 

where is supported on {\x/n — j\ < 1}. Applying this decomposition to both u 
and v, we get 

X n = [ B n f^ V )d V , 

iM% d (10.14) 
/« := f e^e-^-^^^^^^e-r/)^. 
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If j 7^ k, then we choose a G {1, . . . , d} such that |(j — k) a \ ~ \j — k\ and integrate 
the phase e~ in ^- k )'n j n the variable r/ a in K times for arbitrary K E N. Then we get 

/ B n f? k d V = T - k) a ]- K [ [ e UH ^e- m ^ mk ^ 

J a+ p +1 =K JJ (10.15) 

x d«B n dffi(t, V ) (-a B )X(*. Z - v)dtd v . 

Thus denoting 

<P$ ■= ^e-^^iO = n*[-i{z/n - j) a Mx/n - j)<p(z), (10.16) 
we have 

d k r 

„ — n ^ — n 



a=0 a=0 ' 



^oV= ^^(t,x + 2/M(t,x), F ^:=0 (a > 0), 



10.17) 



where in the last term, the summation for (j, k) is constrained to the set {(j — k) a 
\j — k\}, which should be chosen to partition {j ^ k} by a — 1, . . . , d. 
Applying Holder in (£, rj), we have 

d K 



a=0 a=0 

To estimate the L| norm, we recall the standard proof of the (non-endpoint) 
Strichartz estimate from the decay of e . After using Plancherel in (£,77), we 
get 

IK,Jl! L =||| (e < <^-i^ (*,x,y)|l^ (* / )a :,y)> x dydfcft / ) (10.19) 
and the inner product is bounded by using the L p decay of e ltH , 

< it-trii^^^^ikiii^j^^^iiLg, (10.20) 

where 1 < q < 2 and cr = d(l/q — 1/2). Applying the Schwarz inequality to oh/, and 
Hardy-Littlewood-Sobolev to dt'dt, we get 

Rall^^ll^^yJIIirw (10.21) 

for any dual Strichartz exponent (p,q), except for the endpoint, which we do not 
consider for simplicity. By the support property of ttf'" and vj^ (llO.lSp . we have, 
noting that \y\ > n when a > 0, 

IK«(*> x >f)IIW£ <^ a ||min(l,|y/n|- K )«(i,x + y)7j(t,x)|U ?i ^. (10.22) 
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On the other hand, by the Fourier support property of B n , we have 

\\d a qa B n \\ LTL 2 <n a - s \\B n \\ LTti s. (10.23) 
Putting them together, we obtain 

||^n||i| % \\B n \\ L ^H^ S \\(y/ n y Ku ( x + y) v ( x )\\L p t L2 y L% (10.24) 
for any K G N. By using Young's inequality, we have for sufficiently large K, 
\\(y/nr K v(x + yMx)\\« = \\(y/n)- K *{\<p\* * \m{v)\\ L * 



< II (y/n)- Kq \\ Lt ,^ || |</flU/* || (10-25) 

Applying this to the above at each i, then using Holder in t, and summing for all 
n G 2 Z , we get the desired estimate. □ 

10.2. Spatial integration of phase. First we consider the spatially non-resonant 
part Bf. More precisely, we first decompose Bj dyadically by using (|2.6p . 

B]' b ' c := X a (0x b (r/)X C (C)^ = Bf^ x + i?;^ T , (10.26) 

such that |£| ~ a, |^| ~ b and |£| ~ c in the support of Bj ,b,c . The sum over 
a,b,c G 2 Z can be restricted to 

a<&~c, 6< c ~a, c < a ~ 6. (10.27) 

The smooth decomposition into B^' b,c,x and B^ b,c,T will be given in Section [TTJ For 
each Bj' b,c,x , we integrate the phase in 77, by using the identity 

e isn = . ■ V v e lsn . (10.28) 



Then we get terms like 



T' 1 J^ds J e lsn {Bt ■ V V [Z ± ( V )Z ± {^ - v )] + S^fo)^ - r?)} cfy 

= T e iS " {^i[(JZ) ± , Z±] - B^, (JZ)±] + B 2 [Z±, Z±]} ds, 
where f2 = =p H(rj) =F (£ — 77) and £>fc = £>^' c is defined by 



(10.29) 



ja,b,c ■, 
',3 

^ |v,fi| 2 J ' ' 2J ^ |v^| 2 ■ 

Denoting 

M := max(a, b, c), m := min(a, b, c), Z := min(6, c), (10.31) 
we assume that if M -C 1 then 

l|£i Ac ||[^] < ^ 1/2 " 2£ , ll^'l^] < l 1/2 ~ 2£ M~\ (10.32) 
for £ > small, and if M > 1 then 

l|B? 6,c ||[fl3/a+«] S l^ia)- 1 + m- £ , ||B 2 aAc ||^/ 2+£ ] < /^(a)- 1 +m- £ , (10.33) 
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for \e\ small. Note that we have a derivative loss for the first term if (a) C & ~ a 
This is the reason for our new normal form (cf. the discussion in Section [3D- For the 
definition of those norms, see (110.21) and (110.91) . The above estimates will be proved 
in Section [TTJ Then applying Lemma [10.11 with some fixed s = 1 + e > 1 and 3/2, 
we get 

[ e' sH {B 1 [(JZ) ± ,Z ± ]±B 1 [Z ± ,(JZ) ± ]}ds 

a,b,c 

< || r« £ / 2 H?' 6 ' c || [s x + ,] || C/- £ / 2 ^ || ir ^x /2 - £/ 6 1| f/- 1/6 ^ || ^- £ / 4i g 
+ ||^f/(mr/ 2 (a)(6)- 1 (c)- 1 <' c || [B3/2] 



Hi 



(10.34) 



M>X 



X 



\u- £/2 (v)jz\ 



■ l/2-e/6 



\u- 1/6 (v)z\ 



L? /4+e/4 L6> 



where the sums are for a, 6, c satisfying ( 110. 271) . and we used the convention L l l p = 
L p . The norms on B^ b ' c are bounded by the above assumption and interpolation as 
follows. For the first case M < 1, we have 



M<1 



«<1 



(10.35) 



for e' > small, so we can replace the norm with [i5 1+e ] by the interpolation (110. 8p 
for varying e'. For the second case M > 1 we have 

II £ [/^^(a)^)- 1 ^- 1 ^!^/^] < ^[/H^p'" 1 < i, (m36) 

m>i ie2 z 

for |e'| small. Then we can change the norm to [i5 3//2 ] by the interpolation (110. 8p . 
The norm on JZ is bounded by the Sobolev embedding 

\\u- £ i 2 {v)Jz\\ LrLl/ ^ l(i < \\Jz\\ Lrm , 

and the norms on Z are bounded by splitting into \t\ < 1 and |t| > 1, 



(10.37) 



f/- 1/6 z|| LW 4 L6 + ||fr 1/6 (v>z| 



r 3/4+s/4 



LI 



< wu-^zw^e + \\tu- l /*z\\ LrH ^ 6 < \\z\\ snx . 

When restricted for large t > T, the above is bounded by T~ £ / 4 . 
In the same way, we estimate the term with £>2 by 

e isH B 2 [Z ± ,Z ± ]ds 



10.38) 



Hi 



< 



1,6, c i 



(10.39) 



. „ ^ U (l)U (M) l l Q (a) {by 1 {cy 1 B^ || [Bl+£]+[B3/2] 

a,b,c 

X ||f/ _1 Z|| L oo // l||f/ _1/6 Z|| (Ll _ e/2nL 4/3 )tB l,6, 

where the L\ norm is applied to the / frequency component, while the L\ norm is 
applied to the other with frequency max(6, c) ~ M. We use [i5 1+e ] for M < 1 with 
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L] £ ^ 2 L®, and [B 3 ^ 2 ] for M > 1 with L^ 3 L®. Then the norm on £>2 is estimated by 

ll£--A a ' 6 'V+ e ']< ^i 3/2 - 2e 'M- 5 / 6 <i, 

- 1 „ (10.40) 

iiE---^ Ac ii[^ + -]< $>(0<t - x <i, 

M>1 / £ 2 Z 

for e' > small and \e"\ small. Then we can replace the norms with [5 1+e ] and 
[£> 3 / 2 ] respectively, by the interpolation (110. 8p . 

The norms on Z are treated in the same way as for B\. For large t > T, we have 
additional decay factor T~ e l 2 . 

10.3. Time integration of phase. For the temporally non-resonant part Bj ,b,c,T , 
we integrate the phase in s. Then we get the time derivative of the bilinear functions, 
for which we use the equation again. Thus we get terms like 

f e isH Z±] + B 3 [sM±, Z±] + BslZ*, s.A/f ]} cfe 

Jo (11). 41) 

where we put 

£ 3 = B a z f := ^i?; Ac ' T , (10.42) 

with j = 3,4. Recall |£| ~ a, \r}\ ~ 6, |C| ~ c, M := max(a, 6, c), m = min(a, 6, c) 
and Z = min(6, c). We assume that 

11*8*1 [*•] £ ( (M) / M) s l 3/2 ~ s (ay 1 , (10.43) 

for < s < 2. See (110.91) with (110. 2p for the definition of norm. This bound will be 
proved in Section [TT1 

For the first term in (110. 41[) , we use Lemma 110.11 with some fixed s = 1 + e > 1 
for M < 1 and s = 3/2 for M > 1. Thus we get 

< || £f/(/)f7(M) 1 /6 (a)(6) -i (c) -i^ c || [Bl+£]+[i? 3 /2] (1 °- 44) 

Z7 _1/6 (V)Z|| (L l- e /2 nL 4/3 )tL 6||t/ _1 (V)Z|| L oo L 2. 

have 

M<1 1<1 

iiE--- B 3 ,6,c ii [ i,3 /2+e ' ] < 5>wr'- 2 <i, 

M>1 Ze2 z 

for |e'| small. Then by the interpolation (llO.Sp . we can change the norms to [-B 1+£ ] 
and [-B 3//2 ], respectively. If the time interval is restricted on (T, oo), then we get 
additional decay factor T~ £ l 2 from the L\ decay of Z. 



m. 



a,b,c 
X 



For the norm on B3, we have 



(10.45) 
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For the time boundary term, we use Corollary 110.31 with s = 1 and 3/2. Then 



a,b,c 



< || ^t/(/)t/(M) 1 / 3 (a)(/)- 1/2 (M)- 1 ^ c || [i3l]+[B 3 /2] (10.46) 



a,b,c 



X 



lu^HziUeWu-'zw^, 



where we used H 1 C H 1 ^ 2 ' 3 for the lower frequency I term when M > 1. The norm 
on $3 is bounded by 

- 1 , , (10-47) 

iiE-^i^^E^w^ 3/2 <i, 

M>1 ie2 z 

for small, which change into the desired norms by the interpolation. For large 
t > T, we can gain T~ £ by replacing [B 1 ] with [.B 1+E ], and H], with Hx' 1 ~ £ . 

To the remaining terms in (I10.4ip . we apply Lemma 1 10. II with some fixed s = 
1/2 + e > 1/2 and s = 3/2. We obtain 

< || E^(/) 1/2 f/(M)(a)(6)- 1 (c)- 1 ^ c || [Bl/2+£]+[ B3/ 2] (1 °' 48) 

a,b,c 

X ||[/- 1/2 (V)W^|| (I;1 -e/ 2ni2)ti 3||[/- 1 (V)Z|U rX| . 

For the norm on 3%, we estimate 

AKl Kl 

iiE--- B 3 Ac ii [ff a/ 2+e '] < E^)(o^ 2 <i 

for |e'| small, which yields the bound on the desired norm via the interpolation. 

It remains to estimate the above norm on Hz- For the bilinear part of Hz, we 
use (12.201) . splitting into \t\ < 1 and \t\ > 1, 

|| U-^Bjlv*, u ± ]|| i a JT i.a < \M L 2 H ^\Ml?l<>, 

WfU-^B^^W^z < \\tv» 2 

for j — 3, 4. The cubic part of Az is bounded by using 

\\tC j [v ± ,V ± ,V ± ]\\ {L i nL 2 MH l nH l,S )x < WtvWL^HifiWU^vWl^H^Hi^, (10.51) 

for j = 1,2,3,4. Then we apply Sobolev U~ l/2 (H l n F 1,3 ) s C #£> 3 at each t. For 
the quartic term we have 

||Qi(tOII(irW.»). < IMIl« < (tr 2 \H* x , (10.52) 



(10.49) 



(10.50) 
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which gives a better bound than the above. Thus we get the desired bound on J\fz, 
completing the estimate f)10.48p . If we restrict the time interval into (T, oo), then 
we get additional decay T~ £//2 , where the worst decay comes from the bilinear part. 

11. Estimates on the bilinear multipliers with divisors 

In this section we show that we can decompose the bilinear multipliers Bj, j = 3, 4 
such that flT032D . (011331 and (110^31 hold, for all combinations of Z ± , i.e., Bj[Z, Z], 
Bj[Z, Z] and Bj[Z, Z\. We will use elementary geometry in the Fourier space in a 
way similar to [18]. Recall our notation in (12.111) . First we have for |£| > 

|£| - (11-1) 

where H' denotes the radial derivative of H(£). For the higher derivatives, we have 

l Vfc #(OI < |^|=T> (11-2) 
for any k G N, hence by using the Taylor's formula when |£| ~ \rj\, we obtain 

for |£| > \r)\. 



V k H(0 - V k H(rj)\ < ^Z lk _[\ (11.3) 



11.1. Multiplier estimates for ZZ. In this subsection, we deal with the bilinear 
terms of the form 

e UH B J [Z 1 Z]=F- 1 J ^B^rj^-^Z-^Z+^-^drj, (11.4) 

with j = 3,4. The phase Q and its derivatives are given by 
n:=H{£) + H(ri)-H(S-7i), 



vfn = v#(0 - vh(£ - rj), v„n = vfffo) + v#(£ - 77)- 

Recall our notation for the dyadic component 
|f| ~ a, M~6, |C|~c, 

M = max(a, b, c), m = min(a, b, c), Z = min(6, c). 
We also denote in this subsection 



;n.5) 



;n.6) 



or := ic-ei, p--=\c+v\, v ± --=-*;x£xv- (H-7) 

We decompose the (£, 77, £) region (£ = ?7 + C) m to the following five cases smoothly 
and exclusively. Namely, we exclude the earlier cases from the later cases. 

(1) |t7|~|£|^>|C|; Temporally non-resonant. 

(2) a > y/3; Temporally non-resonant. 

(3) |C| ^ I; Spatially non-resonant. 

(4) I?]- 1 - 1 <C M|?y|; Temporally non-resonant. 

(5) Otherwise, spatially non-resonant. 

Now we confirm the desired estimates in each case. 
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11.1.1. \r}\ ~ |£| ^> |C| ; into B T . We put those dyadic pieces with min(6, a) > 8c into 
Bj in ffiTTTT) . In this case we need to show (flTTm Since Q — if (£) + H{rj) — H{Q 
and m = \(\ <C M, we have 

\n\ = Sl>H(M)~M(M). (11.8) 
We choose the (£, () coordinates to use the smallness of ( region. Then 

Vjfi = -V k H(() + V k H(-r]), VfV^fi = -V 1+fe /f (C), (11-9) 



with trivial bounds 



|V c O| < (M>, |V 2 ft| < (m)/m 

\vfn\ < (M), ivjv^ni < (m>/ 



Thus in this dyadic region we obtain 



11.10) 



- — = ™ 3/2 ~ s (M) 1 (11-11) 



_J — r«A c <; _ ,„ 

n 3 



J OO LIS 



M(M) m s (M) 



for j = 3,4 and s — 0, 1, 2, which implies fll0.43[) by interpolation. 

In the remaining cases, we use the coordinates (£, rj) since M ~ \(\ and / ~ |r/|. 

11.1.2. a > \f?>, B T . More precisely, we cut-off the multipliers by 

*[«]:= r(£-c), (H.12) 

for a fixed T G C°°(M 3 ) satisfying T(x) = 1 for \x\ > V3 and T(x) = for \x\ < 3/2. 
Thus we have a > 3/2 in this case, and a < \/3 in the remaining case. By exclusion 
of the previous case, we have M ~ |£| also. Hence we have 

IVJxmI <c^~M" fc , (11.13) 

for any fceN, which are acceptable errors in all of the following cases. 

Now we show (110.431) in this region. Since a > 3/2, we have \rj\ — \(\ > |£|. Hence 
m ~ |£| and 

|Q| = n > H{rj) - H(Q > (M)|e| ~ (M)m, (11.14) 

and from flll.ip . 

iv^ni < ^yH + |v,n| < + (11.15) 

Since a > 3/2, we have j3 < 1/2 and so by the sine theorem (3 ~ aa/6 ~ m/M. 
Thus we obtain 

,_(„)_, . M 2 . . . ._ _ kl . Mm (M)m . Ifi] 

|v " n| ~(M) + <m> ~ (M> ' |v ' n| ~ (m> + m ~ !?• (1116) 

|v,*vf n| = |v 1+t ff(c)| < 

By (111.31) . we have also 

\V 2 Q\ < ^ < H. (11.17) 

^ M 2 ~ M 2 
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Hence we have in this region 



< 



(M) M 3 / 2 m M 3 / 2_s / 3 / 2 - 



(M)m M s (m) (m) 



(a) 



11.18) 



for j = 3, 4, and for s = 0, 1, 2. This yields (110.431) by interpolation. 

11.1.3. |C| > 1, into B x . We put those remaining parts with c > 1 into 5 X . By 
exclusion of the previous cases, we have M ~ |£| > 1 and a < y3. We want to 
show f)10.33p in this region. For the first derivatives, we have from flll.ip . 



ivf^l ~ HCI - leil + |v„n| ~ iiCI - \v\\ + (v)P- 

On the other hand, we have 

Kl 2 = HCI - Ml 2 + 2|dM(i + C • v) ~ IICI - \v\\ 2 + ICIM/? 2 , 
M 2 = IICI - Kll 2 + 2|diei(i - C ■ 9 ~ IICI - l£ll 2 + iciiei« 2 . 



Hence we have 



|v,n|~||ci-HI + (c>/3>|^|. 



For the higher derivatives, we have from (111.31) . 



v;v^n|<(0ICI 



-k 



M 



1-k 



iv^i<ieii^r- fc <iv^i 



\l-k 



11.19) 



11.20) 



11.21) 



11.22) 



Therefore we have in this region 

X[a) B j 



1,1+3/2 

< Z U(a) = l 5 / 2 - s (a)-\ 



ab s 



HS, 



|v„n| 2 

|V„fi| 2 



+ 



Hi 



\v v n\ 



X[ a] Bj 



11.23) 



V,[xg]B. 



<^ U(a) ^^ {ar , 



H* 



for j — 3,4 and s — 0, 1, 2, where X[ a \ = 1 ~ X[a] is excluding the previous case. The 
above estimates imply f|10.33[) by interpolation. 

11.1.4. <C M|?y|, into B T . More precisely, we cut-off the multipliers by 

m := x{(M)\r) x e|/(100M6)), (11.24) 

with x £ Co°(K0 satisfying xil 1 ) = 1 f° r I A 4 1 — 1 an d x(/-0 — for \n\ > 2. Then we 
have 



|V*X[±]I < ((M)/Mb) k , 



11.25) 



which is supported around l^ -1 ! ~ Mb/ (M) for all k > 1. We have included (M) for 
later reuse, but in this subsection we can ignore it because M< 1. More precisely, 
we have by exclusion of the previous cases, 



1 > M ~ |C|, a < v 7 ^. 



11.26) 
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Now we prove (110.43}) in this region. First if M ^ |£|, then Q > H(m) > m and 



iv fc vi 7?) ni < 



M 1 k i V fc+i n i < l£l 



M 



M\7]\ k ' 



for /c > 0, by fll 1 .3j) . Hence we have in this region 

51+3/2 



X[±]X[ a ] B j 



n 



< 



mM(Mb) 



: a~Z 3 / 2 (M/)- 



11.27) 



(11.28) 



for j = 3, 4 and s = 0, 1, 2, where xf a ] * s due ^° exclusion of the second case. Note 
that in this case the smallest divisor (Mb) s is from the cut-off X[±)- The above 
estimate implies (110.431) by interpolation. 

Next we consider the main case M — \(\. Let 



A:= KI + M-ICI- 
Then by elementary geometry, we have 



11.29) 



|^| ~ H(a + /3), A~m(a 2 + /5 2 ) -md^l/l^l) 2 . (11.30) 

Hence the assumption in this case implies that A <C M 2 m. Now we use the small 
but non-zero curvature of H(£) in the radial direction around £ ~ 0, to get a lower 
bound on Q (cf. [M (4.52)]): 

- A = [H(\£\ + \v\)- H(0 - H(rj)} + [H{t -r,)- H(\£\ + M)], 

i^(e-i7)-fr(iei + M)i<<M)A, 

leiMd^i + M) 



11.31) 



Thus we obtain 

For the first derivatives we have 



(£> + (v) 



101 ~ M 2 m. 



VfQ\ < M\r)\ + a « \v v n\ < M\t\ + p < |f |, 



since 



|f |a ~ M/3 ~ m(a +(3) < ml ~ |£||ry|. 
For the higher derivatives we have from fll 1.3j) . 



|vjvf } n| < M~ k , 



lei 



mm* 



Using the volume bound M 2 6 3 as well, we obtain in this case 

bMb 3 / 2 



c 



n 



X[±]X[a] B j 



< 



11 5 



M 2 m{Mby 



l 3/2 (Miy 



11.32) 
11.33) 
11.34) 

11.35) 
11.36) 



for j = 3, 4 and s = 0, 1, 2, which implies (110. 43p by interpolation. 
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11.1.5. Otherwise, into B x . By excluding all of the previous cases, we have now 



1>M~|C|, a<V3, 



> 



M\rj\. 



11.37) 



We want to prove (110.321) in this region. For the first derivatives, we have from 

(mm, 



|V„fi| ~ M||C| - \r}\\ + P, ~ M||C| - + a. 

For the angular part, we have by the sine theorem 

«//3<M/iei, 

since a is away from 2. For the radial part, we have from (111.201) 

Miici-ieii <m\ v \, >mici. 

Thus we obtain 

iv^oi/iv^i < \ v \m. 

Using ( 11 1.35ft for the higher derivatives, we get 



|v„n| 



< 



lei 



< 



i 



|V„fi| 



where we used 

Now we apply the dyadic decomposition to 



< 



< 



\V\ 



M k l3 ~ M^ICHV 



/i G {A; G 2 Z | fc > Mb}. 



Getting a factor /i|^| 1//2 from the L? volume of each dyadic piece, we have 



X[X]X[a] B 3 



h s 



H>Mb 



v,v?"nv,n c c 

X[±]X[a]-°J 



+ 



|V^| 2 



\V v tt\ 



iv^i 4 



X[±]X[a]-Bj 



sE5s«s < 3/! (*«)- 



a/x 



11.38) 



11.39) 



11.40) 



11.41^ 



11.42) 



11.43) 



11.44) 



11.45) 



for j = 3,4 and 1 < s < 2, by using interpolation for each dyadic piece in \i. The 
condition s > 1 is only to have convergence of the sum in /i. As before, xf a ] = 1 — X[a] 
and xf±] — 1 — X[-L] are excluding the previous cases. The above estimates imply 

mm . 
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11.2. Multiplier estimates for ZZ . We consider the bilinear terms of the form 
e itH B j [Z,Z]=F~ 1 J e^Bjirj^-^Z+i^Z+d-n)^, (11.46) 
with j = 3,4, where the phase is given by 

Vfn = VH(0-VH(t- V ), V^ = -Vtf(r/) + Vtf(£-r/). (1L47) 

Recall that |£| ~ a, ~ 6, \(\ ~ c, M = max(a, b, c), m = min(a, b, c) and 
I = min(o, c). By symmetry, we may assume without loss of generality that 

b<c, (11.48) 
in this subsection. We also denote in this subsection 

a:=|f-Cl, P'--=\V-C\, 1--=\Z-V\, tt 1 := -f x f x 77. (11.49) 

We decompose the (£, i], () region into the following four cases exclusively. 

(1) |r/|~|C|^>|^|; Temporally non-resonant. 

(2) a > "\/3; Temporally non-resonant. 

(3) If)- 1 ] -C M\t]\/ (M); Temporally non-resonant. 

(4) Otherwise, spatially non-resonant. 



11.2.1. \r]\ ~ |C| ^> into B T . We put the dyadic pieces with min(o, c) > 4a into 
B T . Then we have M ~ |£| ~ \n\ ^> |£| = m and 

\Q\ > H{rf) ~ M (M), \V k q Vfn\ < (M)M~ k , \V V Q\ < (M)M~ k , (11.50) 

such that in this case 



< < M > M3/2 ™ - ^- s m 
~ M(M) M s (m) M x '' 



11.51) 



for j = 3, 4 and s = 0, 1, 2. By interpolation, this is better than (110. 43j) . 



11.2.2. a > V3, into B T . We use the same cut-off X[a] as defined in ( 111. 12ft . so that 
we have a > 3/2 in this case, and a < v^3 in the remaining case. Then in this region 
we have M — \r}\ ~ |£| ~ |£|, and 

\n\>H(0>M(M), (11.52) 

so that we end up with the same bound as in the previous case. 
In the remaining cases, a < \f?> and \q\ < |C|, hence we have 



(11.53) 
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11.2.3. <C M\rj\/(M), into B T . More precisely, we decompose the multipliers 
by the cut-off X[±] defined in (lll.24p . By exclusion of the previous cases, we have 



now 



M~ |f | ~ |C| > 



For the higher derivatives, we have from ( 111.31) 



|V 1+fc fi| < 



(m) 



m K 



If Kl — If I ^ \v\i then we have 



> (C) - H(0 > (M)m, |V e O| < ^m, |V,0| < (M), 



and so 



-^-X[±]X[al^j 



< 



3/2 



M m 3/2 



11.54) 



11.55) 



11.56) 



11.57) 



for j — 3,4 and s = 0, 1,2, where xf a ] = 1 — X[a] eliminates the previous case. By 
interpolation, it gives (110.431) . 

In the main case |£| — |f | <C \r)\, let A' := |£| + \r]\ — |f |. Since a < 7 and both are 
away from 2, we have 



M 2 

A' ~ m(a 2 + 7 2 ) ~ m 7 2 ~ m{\r) \/\r]\) 2 < J^2 m - 



Then by the same argument as in (lll.3ip . we get 

M 2 



(M) 



m. 



11.58) 



11.59) 



For the first derivatives, we have from (jll.ip . 



ivf } n| - 



M 



(M) 

Combining them with 



we get 



\\C\-\C\\ + (M)a, \v,n\ 



M 

W) 



HCl-hll + W. (11.60) 



m 



a < — 7, 8' ~ 7 
~ M 



IrrH M 
-r-r < 



(M) 



|V^n| < m, |V„ft| < M. 



\2,™3 



Using the volume bound (M/(M)) m as well, we obtain in this region 



< 



(M) Mm 3 / 2 ( (M) \ 



M m 3 / 2 



M 2 (M) \MmJ (M) (M) \Mm J 



(M)V 



11.61) 



11.62) 



11.63) 



for j = 3, 4 and s — 0, 1, 2, which gives (I10.43P by interpolation. 
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11.2.4. Otherwise, into B x . In the remaining case, we have 

M ~ 1^1 ~ ICI > \V\ ~ rn, |^| > Mm/(M). 
Using a < 7 ~ /3' also, we have 



a ol ^ \r\\ m 

r\_< \ 



M 



/3' 7 ~ ICI M' <M) 



< 



|r/- 



»7 



< 7 ~ 



11.64) 



41.65) 



and in particular, 7 > 1 if M > 1. In that case, we have from (111.601) and (111.201) . 



Vf Q\ ~ HCI - ICII + Ma ~ H, |V„fi| > (m). 



If M < 1, we have 



iv^ni < Mm + — ~ — iv„ni. 



M 



M 



Thus in both cases we have 



\vfn\ 
|v^o| 



< 



m(M) 



M(m)' 

Using (111.551) for the higher derivatives, we get 



|vi+*n| 1 



1 



|v„n| 



Thus in this region we obtain 



< 



Vg^O| 

iv^i 



< 



m(M) 



M k (m)\r) A 



41.66) 
41.67) 

41.68) 
41.69) 



|V„ft| 



(M)/im 1/2 M 
~ ^ M(m)fi s JM) 



3/2 



Vj ' v f "v,,si c c 

-X[±]X[a]-Dj 



+ 



|V^| 2 



+ 



vf ) fi(V r? fi) 2 V^ 



c „,c 



IVjjOI 



X[±]X[. 



5, 



Oil J 



(m) 



(11.70) 



IVr,^! 



< 



E 



m{M)iim 1 / 2 M 
M(m)fi 1+S (M) 



.3/2 



(m> 



for j = 3,4 and 1 < s < 2 by interpolation for each dyadic piece of /x, where 
/io := Mm/ (M), the sum for \x is on 2 Z , and xf±] an d Xui are removing the previous 
regions. The above estimates imply (110.321) and (I40.33H . 



11.3. Multiplier estimates for ZZ. For the bilinear terms of the form 

e itH B,[Z,Z] = J e^Bjirj^-^Z-^Z-^-^dr], (11.71) 

we need not decompose, since it is always temporally non- resonant. Indeed we have 

Q = H(0 + H{q) + H(£ - 77) > M(M). (11.72) 
By symmetry, we may assume that \r]\ < \(\. Then 



vf n\ + \Vr,n\ < (M), |v;vf n\ < 



(M) 
M k 



iv; +fc fii < 



41.73) 
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hence we get 



c B 



n 



< <«&,.) < ( 3/-<M)-, 



M(M)b s 



11.74) 



for j = 3, 4 and s = 0, 1, 2, which gives fllQ.43[) via interpolation. 



12. Estimates with phase derivative in cubic terms 

In this section we estimate the phase derivatives in the cubic terms Cj, 1 < 
j < 5. The higher order terms Qj{u) have been already estimated in (19. 8p . by using 
regularity gain. Some parts in the frequency interactions of Cj can be also estimated 
in a similar way, but if we thereby estimate the whole interactions, we can not avoid 
losing regularity by the factor sV^fi. 

On the other hand, the cubic terms have in general a large set of interactions 
which are resonant both in space and time. Hence we cannot always integrate by 
parts as for the bilinear terms. 

Fortunately, it turns out that those two difficulties are disjoint in the Fourier space. 
Namely, when the derivative loss by V^fi becomes really essential, the interaction 
is not temporally resonant, i.e., is even bigger than |Vgf2|. Hence we integrate 
on the phase in s for those parts, and the other parts are bounded in the form V^f2. 
Here we do not need such precise estimates on the resulting multipliers as in the 
bilinear case, since the divisor appear only when |£| ^> 1. 

Before applying J to the cubic terms Cj, we apply the frequency decomposition 
to each v as in (12.61) : 

v =^2 v k, v k :=x k (V)v, ( 12 .1) 
Thus in the Fourier space we get terms like 

£ f sV^e isH Cj[v^v^v%\ds, (12.2) 

where 

Q = H{£) =F #(6) T #(6) T #(6), V c fi = Vtf (0 =F Vif (12.3) 

and one of v has to be replaced with Z if j = 5. We consider the summand with 
k\ > max(/c2, £3), choosing the smaller £2,^3 as integral variables. The other cases 
are treated in the same way. 

If k\ < 1, then all |V.ff(£fc)| are bounded, hence we can bound it in L°°H l by 

\\ tC A v kv v k 2 i v k^\\tfLS/s ^ \\tv\\ L o. L 6\\U' 1 v\\ L o CL 2\\U' 1 v\\ L 2 L 6. (12.4) 

If 1 -c: ki ~ max(/c 2 , ^3), then we have |V^fi| < max(/c 2 , k 3 ). Hence the L°°H l norm 
is bounded by 

\\^X tC A v kv v k^ v k 3 \\\L 2 L^ ^ II^IU-Hi.ell^^lU^^II^^II^Hi.s. (12.5) 

In both cases we get T -1 / 10 for t > T from the V\ decay of U~ 1 v. If j = 5, we just 
replace one of v by Z. 



40 STEPHEN GUSTAFSON, KENJI NAKANISHI, TAI-PENG TSAI 

It remains to estimate the terms with k\ 3> max(l, &2, k^). If the sign in the phase 
is Q = H(£) — H(£i) ± ■ ■ • , then we have 

Vffi = F(£, 6 + 6) • (6 + ^(£, ^ := / V 2 F(e - 9 V )d6, (12.6) 
getting the factor £ 2 + £3. In the region |£| 3> (77}, we have 

Hence the Coifman-Meyer estimate (I2.19P applies to F, and so we can bound the 
Strichartz norm as in the previous case. 

Finally, if Q = + H(£i) ± • • • with &i 3> max(l, fc 2 , £3), we have 

M > kl |V e O| < fcx, (12.8) 

so that we can integrate on e tsn in s. Thus we get terms like 

' V ' Q e™ H {C[v,v,v] + sC[K,v,v]}ds+ ^^„ isH 



n 



n 



-e lsH sC[v,v,v] 







12.9) 



where we omitted indexes j, k* and ±. Moreover if j = 5, one of v (or M v ) should 
be replaced with Z (resp. Nz\ Recall that Af v and Mz are the nonlinearities 
defined in (13.71) and (13.181) . respectively. The multipliers have no singularity thanks 
to the frequency restrictions. Hence using the Coifman-Meyer estimate (I2.19p . we 
can bound the above term in L°°if 1 by 

\\C{v,V,v]\\ L 2 L 6/5 + \\tC[J\f v ,V, v]\\ L 2 L 6/5 +L i/3 L 3/2 + \\tC[v, V , V ) \\ L ™L* 

< \Ml1l4 U ~ 1 v\\ 2 Lool3 + (\\M \\ L 2 L 2 +L 4/3 L 3 + |M|L-L6)ll tl/2f/ ~Mli-LeM 

where the norm on N v is bounded by L 2 ^ 1 ' 6 / 5 + L 4 /3#i,3/2 For J\T Z , it was already 
estimated in Section [HI and the same estimates (18.11) and (I8.3P hold for Af v as well. 
For t > T, we gain T -1 / 5 from the L\ and L\ decay of U~ 1 v (or U~ l Z). 

13. Proof of the initial data result 

Now that we have closed estimates in X for (v, Z) solving the system (I3.16p . we 
can easily finish the proof for Theorem 11.11 More precisely, we have obtained the 
following estimates: For any interval / = (T , T\) C JR., and for any (v, z) G C(I; X) 
satisfying (13. 7p and (13.161) on /, we have 

IKt)-^WIk(*)^<*r 1/8 W)||^ (t) {tell 

h - z \\s(i) < (To)~ r/10 \\v\\ 2 x{I) , 

\\Z-e-^ H Z(T )\\ s(I) < (Tor^dl^l^^ + Ht;!!^^}, (13.1) 
\\Z-e-^- T ^ H Z(T Q )\\ x{I) 

~ (^0) 6 {\\ v \\\xc]S)(i) + IMIpcnsx/) + Il^ll(xn5)(/)} > 

for some small e > 0, where the highest order 6 comes from the partial integration 
replacing one of v in Cj by the quartic term Q\{u). 
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There are at least two ways to procede: (1) construct global solutions with the 
desired properties by an iteration argument, or (2) use those estimates as a priori 
bounds for the global solutions already known in H . 

Since the equations for the iteration sequence are a bit complicated due to the 
partial integration and the nonlinear transforms, we employ the latter (2) approach. 
Then we just need to ensure that the global solutions from our initial data stay 
in X so that the a priori bounds can be rigorously applied. Thanks to the global 
wellposedness [3] in if) G C(M; 1 + i/ 1 (lR 3 )) for (II. ip . it suffices to prove that only for 
a subset of initial data which is dense with respect to the H 1 norm for u. 

It is quite standard to see that for any u(0) G H 2 n (x) H 1 , the unique global 
solution of (13.11) stays in the same space: since we already have it in L°°H l , we just 
apply the fixed point argument to the following estimates on (0, T) 

ll M IU°°//2 < \\Nu\\l2hw/ b +lihi < {T 1/2 + T)\\u\\ L ^ H 2(l + \\u\\ L ^ H i) 2 , 

\\{x + 2itV)u\\ H i < \\{x + 2itV)M u \\ L 2 H i,e/ 5+L x H i (13.2) 

< (T 1/2 + T 2 )(\\xu\\ L ^ H i + \\u\\ L vo H2 )(l + \\u\Ilcvhi) 2 , 

where we denote the nonlinearity in the equation (jl.fip by 

N u ■= u 2 + 2\u\ 2 + \u\ 2 u. (13.3) 

The above estimates follow from Strichartz applied to (11.61) . and the identity x + 
2iiV = e ttA xe~ ttA . Thus we get the solution locally in this space, and then extend 
it globally by Gronwall. 

To those solutions u G H 2 n (x)~ l H l satisfying (11.121) with v — u\ + iUui and 
Z — v + b(u), we can apply the above a priori estimates on any interval, thereby 
getting global bounds 

\Mxns + H^IUns < IMO)Hx(o) < (13.4) 

provided 5 is sufficiently small. By a density argument, all solutions satisfying our 
initial assumption can be obtained as limits of such solutions. Hence the above 
global bounds are transferred to them, and we get the desired scattering for them, 
using the identities 



\e ltH Z{t) - e iToH Z(T )||^i = \\Z(t) - e~ l{t ~ To)H Z{T { 



{e itH Z{t) -j ToH Z{T )}\\ Hl = \\j{Z(t) -e-^~ To)H Z(T )} 



//! ' (13.5) 



\h 1 



Denote the map W : v (0) ^ v + thereby obtained (the inverse wave operator) 
locally around in (x)~ H 1 . The continuity and the invertibility of W are also 
proved in a standard way. The above estimates are adapted directly to differences 
of solutions, since after writing the system (13.161) for (v,Z), we did not use any 
genuinely nonlinear structure, but only multi-linear estimates. Let (v^\Z^') be 
small solutions of fl3~16l) obtained above, and let $ ] = e - i(t ~ T)H Z (j) (T), for j = 1,2 
and T > 0. Then we have 

\\zv - ZP> - - $>)\\ xna 

2 (13 6) 

< ^(ll^ll^ + II^IUn^CII^-^IUns+ll^-^lkns), 
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(13.8) 



and also 

- (Z« - Z^)\\ xns < \\v U) \\xns\\v W -v^\\ XnS . (13.7) 

i=i 

Hence by using the smallness assumption, we obtain 

\\zu - Z^\\ xns ~ - v V>\\ Xn8 ~ \\f® - f?\\xns 

~ ||Z«(T) - z( 2 )(T)|| x(r) ~ || w W(r) - ^ 2 )(T)|U(T). 

Since Pv> (j) (0)) = lim^ e itH Z^\t) in (x^H 1 , the bi-Lipschitz property of W 
follows from the above equivalence. 

Concerning the invertibility, or construction of the wave operator, we can argue as 
for Theorem ll.2[ namely consider an approximating sequence of solutions (v T , Z T ) 
with Z T (T) = e~ %TH v + for any given small v + G (x)~ H 1 , and let T — > oo. Then 
the above a priori bound implies weak convergence of some subsequence in X D S, 
whose limit has the desired scattering property. □ 

Appendix A. Correction for the 2D asymptotic profile in [T8] 

The asymptotic profile given in [18] for the two dimensional case was incorrect. 
The equation for z := U~ l [u\ + (2 — A)" 1 ^ 2 ] + ra 2 was given by0 

iz-Hz = 2u\ + 4iV(2 - A)- 1 U- 1 V(u 1 Vu 2 ) + 0(u 3 ) 

= 2[U(z + z)/2} 2 + tH- 1 V[U{z + z)-V{z-z)}+0{u 3 ), (A " l) 

so the correction term z' coming from z~z should be given by 

I e ^(»-*)||^0|2 +2iImff -i v ^^o. dSj ( A 2 ) 

where z° := e~ ltH ip is the given free solution. The second term was missing in [T8j . 
Nevertheless the proof was correct because it was carried out assuming that all the 
bilinear terms with z°z° had been subtracted, and those correction terms did not 
really affect the estimates. Both terms can give non-L 2 contribution around |£| ~ 0. 
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